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the  molecules  in  a  unit  cell  (or  larger  piece  of  crystal)  in  the  multipole  field  of  yet  fur¬ 
ther  out  surrounding  molecules.  The  SCF  wave  function  is  localized  and  the  localized  orbitals 
(occupied  and  virtual)  in  the  region  of  the  bond  being  broken  are  included  explicitly  in  the 
MRD-CI  calculations.  This  method  is  completely  general  and  applicable  to  any  molecule  in  any 
kind  of  a  crystal  or  other  solid  environment.  This  development  has  led  to  an  important 
breakthrough  which  will  lead  to  crucial  understanding  of  the  initiation  of  detonation  and  the 
subsequent  processes  leading  to  detonation.  Results  will  presented  on  the  CH,  -  NO-  decom¬ 
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including  the  multipole  effects  of  farther  out  molecules  to  include  long  range  effects  also. 

We  then  meshed  this  multipole  procedure  back  into  the  MRD-CI  programs  to  enable  us  to  include 
multipole  effects  when  breaking  a  chemical  bond  in  a  crystal.  We  also  derived  and  implemented 
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We  carried  out  ab-initio  crystal  orbital  calculations  on  several  unit  cells  of 
nitromethane  to  verify  that  our  SCF  method  in  the  field  of  multipoles  described  above  in  Part 
III  did  correctly  describe  a  crystal  of  nitromethane. 
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CONCISE  SUMMARY 


I.  Program  Enhancements  and  New  Program  Developments  on  the  CRAY 

Supercomputer 

This  past  year  we  have  made  a  significant  breakthrough.  We  developed 
and  implemented  and  used  successfully  the  strategy  for  ab-initio  MRD-CI 
(mul tireference  double  excitation  -  configuration  interaction)  calculations 
for  breaking  a  chemical  bond  in  a  molecule  in  a  crystal  or  other  solid 
environment.  In  this  procedure  the  SCF  is  solved  explicitly  for  the 
molecules  in  a  unit  cell  (or  larger  piece  of  crystal)  in  the  multipole  field 
of  yet  further  out  surrounding  molecules.  The  SCF  wave  function  is 
localized  and  the  localized  orbitals  (occupied  and  virtual)  in  the  region  of 
the  bond  being  broken  are  included  explicitly  in  the  MRD-CI  calculations. 
This  method  will  be  detailed  in  section  III.  This  method  is  completely 
general  and  the  results  lead  to  an  understanding  of  fractoemission  and  of 
the  initiation  of  detonation  and  the  subsequent  processes  leading  to 
detonation. 

We  continue  to  serve  on  the  NSF  San  Diego  Supercomputer  Center  (SDSC) 
computer  time  allocation  committee  on  their  CRAY  XMP  4/8.  We  also  receive 
computer  time  grants  from  SDSC  which  we  use  for  the  bulk  of  our  quantum 
chemical  calculations  on  this  ONR  research. 

II.  MRD-CI  Calculations  for  Cationic  Polymerization  of  Energetic  Oxetanes 

Our  major  emphasis  this  past  year  has  been  to  carry  out  in-depth 
detailed  ab-initio  MRD-CI  (mul ti reference  double  excitation  -  configuration 
interaction)  calculations  on  the  propagation  step  of  cationic  polymerization 
of  prototype  substituted  energetic  oxetanes. 

Cationic  polymerization  consists  essentially  of  two  major  steps: 
initiation  and  then  propagation.  There  is  considerable  Navy  interest  in 
energetic  polymers  made  by  cationic  polymerization  of  oxetanes  substituted 
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or  di substituted  by  exotic  energetic  substituents  such  as  azido, 
azidomethyl,  nitrato,  nitraminomethyl ,  etc.  as  well  as  fluoro  and  nitro 
groups.  The  initiation  step  (which  is  crucial  for  cationic  polymerization 
to  take  place)  is  governed  by  the  propensity  of  the  substituted  oxetane  to 
undergo  protonation.  Our  previous  ab-intio  quantum  chemical  SCF 
calculations  on  the  energetic  oxetane  monomers  and  electrostatic  molecular 
potential  contour  (EMPC)  maps  we  generated  from  these  electronic  wave 
functions,  which  predict  the  order  of  protonation  and  hence  initiation,  were 
able  to  predict  correctly  the  propensity  of  the  energetic  substituted 
oxetane  monomers  to  undergo  polymerization  even  prior  to  the  synthesis  of 
the  monomers. 

A.  Ab-Initio  MRD/CI  Calculations  for  the  Propagation  Step 
1.  Discussion 

As  was  suggested  to  us  by  several  different  experimentalists  in 
cationic  polymerization  (primarily  Gerry  Manser)  the  mechanism  seems  to  be 
attack  of  protonated  oxetanes  on  oxetanes  (or  vice  versa)  with  concomitant 
ring  opening  of  the  protonated  oxetane  according  to  the  following  general 
scheme 
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We  have  carried  out  this  past  year  and  are  continuing  to  carry  out  ab- 
initio  MRD-CI  calculations  on  the  subsequent  propagation  step  of  oxetane  (or 
an  energetic  substituted  oxetane)  reacting  with  protonated  oxetane  (or  a 
protonated  energetic  substituted  oxetane). 

MRD-CI  calculations  along  the  potential  energy  surfaces  have  been 
carried  out  for  a  very  large  number  of  geometry  variations  for  the  angles 
between  the  planes  of  the  substituted  oxetane  and  protonated  substituted 
oxetane  rings  (which  can  be  different  in  each  direction  in  the  case  of 
substituted  rings),  the  inter-ring  distance  ( 0 ^ a)  (w^ere  the  A  ring  is 

the  protonated  ring  and  the  B  ring  is  the  non-protonated  ring),  the  angle  of 
opening  the  C^-0^  ring  and  the  orientation  of  the  H  atoms  on  C4^. 

The  preferred  direction  of  attack  appears  to  be  the  reaction  of  the 
oxygen  of  the  unprotonated  oxetane  ring  (which  we  call  O^g)  with  the  a 
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carbon  (which  we  call  C^)  of  the  protonated  substituted  oxetane  ring  along 
the  C^-0^  bond  direction  with  conccmmitant  pulling  back  (inversion  of  the 
H  atoms  on  C^)  and  opening  of  the  C^-0^  bond  in  the  protonated  oxetane 
ring  and  formation  of  an  bond. 

The  ab-initio  MRD-CI  calculations  on  the  propagation  step  of  the 
protonated  oxetane  ring  opening  in  the  course  of  interaction  with  oxetane 
were  carried  out  based  on  localized  orbitals  on  the  pertinent  regions 
involved  in  the  reaction. 

These  MRD-CI  calculations  have  enabled  us  to  map  out  the  reaction 
coordinates  of  the  propagation  step  reaction  of  oxetane  (or  an  energetic 
substituted  oxetane)  reacting  with  protonated  oxetane  (or  with  a  protonated 
energetic  substituted  oxetane),  to  identify  the  transition  state  of  the 
propagation  step  and  to  identify  when  the  bond  in  the  protonated 

ring  will  start  to  open  as  a  function  of  inter-ring  distance  and  angle  for 
each  different  pair  of  substituted  reactants. 

2.  Detailed  Results 

This  year  we  first  carried  out  such  detailed  MRD-CI  calculations 

for  the  prototype  systems  OXET  +  QXETH+,  OXET  +  FN0XH+,  FNOX  +  0XETH+,  FNOX 

+  FN0XH+,  to  gain  the  understanding  of  the  basic  mechanism  of  the 
propagation  step  involving  energetic  substituted  oxetanes  without  the 
additional  complication  of  floppyside  chain  groups.  We  then  carried  out 

such  detailed  MRD-CI  calculations  for  the  systems  AMMO  +  0XETH+,  OXET  + 

AMM0H+  and  AMMO  +  AMM0H+.  We  are  currently  carrying  out  similar 

calculations  involving  BAMO  and  BAM0H+  and  will  carry  such  calculations  out 
for  other  energetic  substituted  oxetanes.  These  results  will  enable  us  to 
understand  and  to  predict  copolymerization  propensities. 

B.  Ab-Initio  MRD-CI  Calculations  of  the  Protonation  of  Oxetane 

Oxetane  +  H+  is  not  the  lowest  energy  state  of  separated  fragments 
of  protonated  oxetane  at  the  dissociation  asymptote.  The  lowest  energy 

state  is  oxetane+  +  H  since  the  IP  of  oxetane  is  lower  than  that  of  H. 

Thus,  no  single  determinant  SCF  calculation  for  protonation/deprotonation 
can  describe  the  system  properly.  Our  MRD-CI  results  indicated  that  the 

lowest  ground  *A^  state  at  equilibrium  dissociated  to  oxetane+  (^Aj)  +  H. 

The  potential  surface  arising  from  oxetane  (X  *A^)  +  H+  was  repulsive. 

There  were  also  a  wealth  of  other  states  arising  at  the  dissociation 
asymptote  from  higher  states. 


\  V  V« 
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III.  Ab-Initio  MRD-CI  Calculations  for  Breaking  a  Chemical  Bond  in  a 

Molecule  in  a  Crystal  or  Other  Solid  Environment 

The  challenge  arose  to  extend  our  MRD-CI  (mul Preference  double 
excitation  -  configuration  interaction)  technique  based  on  localized/local 
orbitals  to  the  breaking  of  a  chemical  bond  in  a  molecule  in  a  crystal  (or 
other  solid  environment).  This  past  year  we  have  derived,  implemented,  and 
used  successfully  a  procedure  for  doing  this.  We  made  the  first 
presentation  of  results  using  this  method,  spring  1987,  at  the  ONR  WorkshoD 
on  Dynamic  Deformation,  Fracture  and  Transient  Combustion  of  Energetic 
Compounds . 

This  development  has  led  to  an  important  breakthrough  which  will  lead 
to  crucial  understanding  of  fractoemi ssion  and  of  the  initiation  of 
detonation  and  the  subsequent  processes  leading  to  detonation.  Our  method 
is  completely  general  and  applicable  to  any  molecule  in  any  kind  of  a 
crystal  or  other  solid  environment.  The  crystal  can  have  defects, 
deformations,  dislocations,  impurities,  dopants,  edges  and  surface 
boundaries,  etc. 

Results  will  presented  on  the  CH^  -  NO^  decomposition  of  nitromethane 
in  nitromethane  crystal.  This  system  is  the  prototype  of  >C-N02 
dissociation. 

IV.  POLY-CRYST 

POLY-CRYST  is  the  program  we  previously  derived  and  wrote  for  ab-initio 
calculations  on  crystals  and  polymers  using  the  translational  symmetry  in  a 
crystal  and  the  translational-rotational  symmetry  in  a  polymer.  Commensurate 
with  the  ONR  priorities  expressed  to  us  by  our  ONR  Contract  Monitor,  we 
devoted  only  minimal  but  still  scientifically  significant  effort  to  further 
development  and  testing  of  the  POLY-CRYST  program.  As  options  we  had  already 
included  in  POLY-CRYST  our  own  ab-initio  MODPOT  (ab-initio  effective  core 
model  potentials)  and  VRDDO  (charge  conserving  integral  prescreening 
evaluation)  options.  It  is  these  features,  particularly  VRDDO,  which  will 
enable  POLY-CRYST  to  handle  molecular  crystals  of  large  molecules  and  with 
large  numbers  of  large  molecules  per  unit  cell.  This  year  we  derived  and 
incorporated  into  the  POLY-CRYST  including  the  multipole  effects  of  farther 
out  molecules  to  include  long  range  effects  also.  We  then  meshed  this 
multipole  procedure  back  into  the  MRD-CI  programs  to  enable  us  to  include 
multipole  effects  when  breaking  a  chemical  bond  in  a  crystal.  We  also  ran 
some  tests  on  POLY-CRYST  on  integral  thresholds  and  numbers  of  unit  cells 
necessary  for  convergence.  These  preliminary  tests  identified  necessary 
criteria. 

Toward  this  convergence  criteria  goal,  we  also  derived  and  implemented  a 
procedure  for  calculating  the  charge  imbalance  caused  by  various  integral 
thresholds  to  give  a  precise  measure  of  the  effect  on  the  crystal  orbital 
calculation  of  dropping  integrals  of  various  sizes.  The  POLY-CRYST  program 
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I.  Program  Enhancements  and  New  Program  Developments  on  the  CRAY 
Supercomputer . 

This  past  year  we  have  made  a  significant  breakthrough.  We  developed 
and  implemented  and  used  successfully  the  strategy  for  ab-initio  MRD-CI 
(multireference  double  excitation  -  configuration  interaction)  calculations 
for  breaking  a  chemical  bond  in  a  molecule  in  a  crystal  or  other  solid 
environment.  In  this  procedure  the  SCF  is  solved  explicitly  for  the 
molecules  in  a  unit  cell  (or  larger  piece  of  crystal)  in  the  multipole  field 
of  yet  further  out  surrounding  molecules.  The  SCF  wave  function  is 
localized  and  the  localized  orbitals  (occupied  and  virtual)  in  the  region  of 
the  bond  being  broken  are  included  explicitly  in  the  MRD-CI. 

This  method  will  be  detailed  in  Section  III. 

This  method  is  completely  general  and  the  results  lead  to  an 
understanding  of  fractoemission  and  the  initiation  of  detonation  and  the 
subsequent  processes  leading  to  detonation. 

We  continue  to  serve  on  the  NSF  San  Diego  Supercomputer  Center  (SDSC) 
computer  time  allocation  committee  on  their  CRAY  XMP  4/8. 

Dr.  Kaufman  served  at  the  December  1986  Allocation  committee  meetings. 
Dr.  Hariharan  served  at  the  March  and  June  1987  meetings  and  Dr.  Koski  at 
the  September  1987  meeting. 

We  also  receive  computer  time  grants  from  SDSC  which  we  use  for  the 
bulk  of  our  quantum  chemical  calculations  on  this  ONR  research. 


II.  MRD-CI  Calculations  for  Cationic  Polymerization  of  Energetic  Oxetanes 


Our  major  emphasis  this  past  year  has  been  to  carry  out  in-depth 
detailed  ab-initio  MRD-CI  (mul ti-reference  double  excitation-configuration 
interaction)  calculations  on  the  propagation  step  of  cationic  polymerization 
of  prototype  substitued  energetic  oxetanes.  Cationic  polymerization 
consists  essentially  of  two  major  steps:  initiation  and  then  propagation. 
There  is  considerable  Navy  interest  in  energetic  polymers  made  by  cationic 
polymerization  of  oxetanes  substituted  or  di substi tuted  by  exotic  energetic 
substituents  such  as  azido,  azidomethyl,  nitrato,  ni traminomethyl ,  etc.  as 
well  as  fluoro  and  nitro  groups.  The  initiation  step  (which  is  crucial  for 
cationic  polymerization  to  take  place)  is  governed  by  the  propensity  of  the 
substituted  oxetane  to  undergo  protonation.  Our  previous  ab-intio  quantum 
chemical  SCF  calculations  on  the  energetic  oxetane  monomers  and 
electrostatic  molecular  potential  contour  (EMPC)  maps  we  generated  from 
these  electronic  wave  functions  which  predict  the  order  of  protonation  and 
hence  initiation,  were  able  to  predict  correctly  the  propensity  of  the 
energetic  substituted  oxetane  monomers  to  undergo  polymerization  even  prior 
to  the  synthesis  of  the  monomers. 

A.  Ab-Initio  MRD/CI  Calculations  for  the  Propagation  Step 

1.  Discussion  of  Calculation  Procedure  and  Pathways  of  Attack 


As  was  suggested  to  us  by  several  different  experimentalists 
in  cationic  polymerization  (primarily  Gerry  Manser)  the  mechanism  seems  to 
be  attack  of  protonated  oxetanes  on  oxetanes  (or  vice  versa)  with 
concomitant  ring  opening  of  the  protonated  oxetane  according  to  the 
following  general  scheme 
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We  have  carried  out  this  past  year  and  are  continuing  to  carry  out  ab- 
initio  MRD-CI  calculations  on  the  subsequent  propagation  step  of  oxetane  (or 
an  energetic  substituted  oxetane)  reacting  with  protonated  oxetane  (or  a 
protonated  energetic  substituted  oxetane). 


In  order  to  understand  and  to  be  able  to  predict  copolymerization 
propensities  of  various  energetic  substituted  oxetanes  it  is  necessary  to 
trace  the  reaction  pathways  of  the  propagation  step  in  cationic 
polymerization . 


Although  geometries  of  reactants  and  products  may  generally  be  obtained 
experimentally  using  a  wide  varietyof  spectroscopic  methods  these  same 
techniques  provide  little  information  about  reaction  pathways.  Thus,  the 
only  way  that  the  pathways  of  reactions  and  geometries  of  reaction 
transition  structures  may  be  obtained  is  from  the  quantum  chemical 


\/ 


calculations  theory.  Such  quantum  chemical  theory  can  be  used  to  examine 
any  arrangement  of  nuclei. 

These  energetic  oxetanes  are  large  molecules  for  MRD-CI  calculations 
and  the  systems  of  energetic  oxetanes  and  protonated  oxetanes  are  even 
larger  and  thus  beyond  the  size  in  which  MRD-CI  calculations  can  be  carried 
out  in  the  cpu  memory  and  disc  storage  of  current  CRAY  XMP  supercomputers . 
Thus  we  had  derived,  implemented  and  tested  a  new  computational  strategy  for 
MRD-CI  calculations  for  intermolecular  reactions  and  for  molecular 
decompositions  based  on  localized  orbitals.  (The  strategy  is  described  in 
more  detail  later  in  this  section). 


MRD-CI  calculations  along  the  potential  energy  surfaces  have  been 
carried  out  for  a  ve^y  large  number  (at  least  25  separate  MRD-CI  calculated 
points  at  different  geometries  are  necessary  for  each  set  of  reacting 
partners)  of  a  angles  between  the  planes  of  the  substituted  oxetane  and 
protonated  substituted  oxetane  rings  (which  can  be  different  in  each 
direction  in  the  case  of  substituted  rings),  the  inter-ring  distances  (O^g- 

C^.)  (where  the  A  ring  is  the  protonated  ring  and  the  8  ring  is  the  non- 


protonated  ring),  angles  6  of  opening  the  C^.-O^ 


bond  in  ring  A  and  the 


orientation  of  the  H  atoms  on  C^A  as  a  function  of  the  inter-ring  distance. 


For  the  MRD-CI  calculations  on  oxetanes  and  protonated  oxetanes 
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we  considered  the  localized  bonds  in  the 


bond,  the  ^ja~^4A 


bond,  the 


^2aXa  dond»  t'he  OiA-H+  bond,  the  lone  pairs  on  °1A  and  the  bonds 

connecting  hydrogens  to  C2^  and  C4^,  the  O^g-C^g  bond,  the  ^ i B~^2 B  donc1,  the 

bonds  connecting  hydrogens  to  and  ^45  and  the  inter-ring  bond. 

This  choice  of  localized  orbitals  has  the  great  advantage  since  energetic 
oxetanes  are  substituted  in  the  3  position  that  it  preserves  the  similarity 
in  the  MRD-CI  among  all  the  energetic  substituted  oxetanes  and  protonated 
oxetanes  and  provides  a  sound  basis  for  comparison. 


The  preferred  direction  of  attack  appears  to  be  the  reaction  of  the 
oxygen  (which  we  call  O^g)  of  the  unprotonated  oxetane  ring  on  the  a  carbon 
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(which  we  call  C^)  of  the  protonated  substituted  oxetane  ring  along  the 
C4A-O1A  bond  direction  with  concomitant  pulling  back  (inversion  of  the  H 
atoms  on  C^)  and  opening  of  the  bond  in  the  protonated  oxetane 
ring,  similar  to  an  S^2  reaction  mechanism. 


\  / 
C 


.  /  . 

,c  "B"  ring 

\  / 

0 

cx>r$  ;  c  "A" 

1 

H 


ring 


The  angle  a  between  the  two  rings  is  determined  for  R(0jg_C4^)  ■  2.6  - 
2.9  -  3.4  bohrs  and  is  used  for  all  other  geometries. 


Two  bonds  are  essential. 


C4A-O1A  and  are  essential  to  describe 


the  reaction  pathway.  The  bond  inside  the  protonated  oxetane  ring  ( A~0 ^ A ) 

varies  from  R  -  2.8  to  R  *  4.9  bohrs  which  correspond  to  the  fully  closed 
and  fully  open  ring  of  the  protonated  oxetane.  This  bond  is  described  by 

O  O 


the  parameter  6  which  varies  from  0  (fully  closed  ring)  to  19  (fully  open 
ring)  and  corresponds  to  the  degree  of  openness  of  the  ring. 


The  01B-C4A  bond  is  changed  from  R ( 0 ^ g-C^ a )  ^rom  R  *  2.1  bohrs  to  R  * 
10.0  bohrs. 


Positions  of  the  proton  H+  and  hydrogens  connected  to  C4A  atom  are  the 

most  affected  by  opening  the  ring  and  their  positions  were  determined  for 

the  prototype  system  OXET  +  FN0XH+.  These  proton  and  hydrogen  positions 
were  then  used  for  subsequent  studies  of  propagations  reactions  involving 
other  protonated  energetic  oxetanes. 

We  had  previously  shown  that  ab-initio  M0DP0T/VRDD0  MRD-CI  calculations 
for  oxetanes  and  protonated  oxetanes  gave  energy  differences  and  MRD-CI 
coefficients  very  close  to  those  from  much  larger  basis  set  all-electron 
MRD-CI  calculations. 


Ab-initio  M0DP0T/VRDD0  MRD-CI  calculations  have  been  carried  out  for 
each  point  of  the  potential  surface  of  oxetanes  reacting  with  protonated 
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oxetanes  in  the  propagation  step  of  cationic  polymerization.  Because  of  the 
size  of  the  intermolecular  complex  molecular  orbitals  selected  from 
localized  space  are  used  in  the  MRD-CI  calculations.  The  geometries  studied 
include  the  most  sensitive  part  of  the  complex  in  the  the  MRD-CI  procedure. 
Ten  of  the  most  important  main  reference  configurations  have  been  used  in 
MRD-CI  treatment,  and  the  same  set  of  main  reference  configurations  have 
been  kept  through  whole  potential  surface.  All  single  and  double 
excitations  were  allowed  relative  to  these  main  configurations.  The 
energies  of  each  of  the  thousands  of  contributing  configurations  is 
estimated  by  a  perturbation  procedure;  a  threshold  is  set  for  which 
contributions  wi 1  be  included  explicitly  in  the  MRD-CI,  in  the  following 
tables,  this  MRD-CI  energy  is  designated  Cl.  Then  the  energies  of  all  of 
the  other  configurations  generated  but  not  included  explicitly  in  the  MRD-CI 
are  extrapolated  and  added  back  in,  this  energy  is  designated  EX.  Finally  a 
Davidson  type  correction  (which  has  been  shown  to  be  a  good  correction)  for 
size  extensivity  is  added  in. 

ref 

E( ful 1  Cl  estimate)  =  E( EX)  +  (1  -  Z  Cp2)  [E( EX)  -  E ( Ref) ] 

P 

and  the  summation  is  over  all  reference  species.  The  use  of 
multiconf igurational  scheme  is  to  assure  avoiding  of  possible  artifacts. 

Our  MRD-CI  results  support  the  suggestions  of  Gerry  Manser  as  to  the 
mechanism  of  the  propagation  step  in  cationic  polymerization  of  oxetanes. 

We  discussed  this  with  Gerry  and  he  was  quite  gratified  that  our  theoretical 
results  were  in  accord  with  his  hypothesis. 

These  MRD-CI  calculations  have  enabled  us  to  map  out  the  reaction 
coordinates  of  the  propagation  step  reaction  of  oxetane  (or  an  energetic 
substituted  oxetane)  reacting  with  protonated  oxetane  (or  with  a  protonated 
energetic  substituted  oxetane),  to  identify  the  transition  state  of  the 
propagation  step  and  to  identify  when  the  C^-0^  bond  in  the  protonated 

oxetane  will  start  to  open  as  a  function  of  inter-ring  distance  and  angle 
for  each  different  pair  of  substituted  reactants. 

By  comparing  these  results  for  different  pairs  of  reacting  substituted 
oxetanes  and  protonated  substituted  oxetanes  we  shall  be  able  to  gain 
insight  into  preferred  copolymer  candidates  and  relative  reactivity  ratios. 

As  we  mentioned  above,  the  substituted  oxetanes  are  large  molecular 
systems  and  their  interactions  with  protonated  substituted  oxetanes  lead  to 
even  larger  systems.  These  are  larger  molecular  systems  than  have  ever  been 
calculated  with  MRD-CI  methods  and  also  exceed  the  amount  of  data  that  can 
be  handled  even  in  the  current  CRAY-XMP  series  with  the  available  core  and 
disc  space. 

We  had  previously  developed  and  validated  a  new  MRD-CI  approach  based 
on  localized  orbitals  in  the  reaction/interaction  region  with  the  remainder 
of  the  non-participating  localized  occupied  molecular  orbitals  being  folded 
into  an  effective  Cl  Hamiltonian.  We  had  shown  by  test  examples  that  the 


MRO-CI  based  on  localized  orbitals  give  a  potential  energy  surface  for 
molecular  decomposition  essentially  parallel  to  that  using  the  entire 
valence  space  MRD-CI.  These  MRO-CI  calculations  for  the  reaction  of 
substituted  protonated  oxetanes  with  substituted  oxetanes  are  a 
computationally  and  labor  intensive  project.  For  each  different  inter-  and 
intra-molecular  geometry  point,  first  the  SCF  calculation  must  be  run,  then 
the  resulting  SCF  canonical  delocalized  molecular  orbitals  must  be 
localized.  In  addition  to  determining  the  centroids  of  the  localized  bonds 
compared  to  the  bonds  and  atoms  involved  in  the  interaction/reaction,  this 
past  year  we  supplemented  this  with  additional  procedures  to  ensure  that  all 
pertinent  bonds,  lone  pairs,  vacant  orbitals  etc.  would  be  included  in  the 
localized  orbital  MRD-CI.  After  localization  a  small  single  reference  Cl 
must  be  carried  out  to  determine  the  major  reference  configurations  to 
include  in  the  subsequent  MRD-CI.  A  great  advantage  in  our  carrying  out  the 
ab-initio  MRD-CI  calculations  based  on  the  important  localized  orbitals  in 
the  interaction/reaction  region  is  the  reasonable  similarity  of  types  of 
major  reference  conf igurations  for  the  variously  energetic  substituted 
oxetanes  and  energetic  substituted  protonated  oxetanes. 

By  first  carrying  out  detailed  MRD-CI  calculations  on  the  partners 
(oxetanes  and  protonated  oxetanes)  of  oxetane  itself,  with  the  prototype 
energetic  oxetane,  FNOX  (3-fluoro-3-ni trooxetane)  we  have  been  able  to 
identify  the  most  pertinent  portions  of  the  interaction  surface  and  their 
electronic  characteristics .  This  behavior  was  verified  by  our  similar 
calculations  involving  protonated  with  non-protonated  partners  of  oxetane 
and  AMMO  (3-azidomethyl-3-methyloxetane)  We  can  now  focus  on  these  regions 
for  our  calculations  involving  other  even  larger  energetic  substituted 
oxetanes  and  protonated  substituted  oxetanes  to  identify  copolymerizaton 
propensities  and  optimal  copolymer  candidates. 
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2.  Detailed  Results  for  Prototype  Examples  are  presented  in  the 
fol lowing  pages. 

(1).  OXET  +  0XETH+,  OXET  +  FN0XH+,  FNOX  +  0XETH+,  FNOX  + 
FN0XH+ 

(a).  Energies 

1'.  Oxetane  (OXET)  +  protonated  oxetane 
(0XETH+) 

a'.  Determination  of  The  Pathway 

The  addition  of  oxetane  (OXET)  to 

protonated  oxetane  (0XET+)  has  the  features  of  an  S^2  reaction  as  we 

reported  preliminarily  last  year.  The  approach  presented  utilizes  these 
features  to  the  greatest  extent.  Geometry  of  the  reacting  system:  The 

reaction  of  OXET  with  0XETH+  is  assumed  to  be  an  S..2  reaction  with  molecule 

N 

B  (OXET)  attacking  molecule  A  (0XETH+)  along  the  0ib-C4/\-0j/\  line-  It  is 

found  that  head  on  geometry  for  approaching  oxetane  is  energetically 
preferred,  thus  the  Cgg  is  also  kept  on  the  0ib_^4A-01A  1ine- 

The  geometry  for  the  oxetane-protonated  oxetane  complex  is  presented 
in  Figure  II-l,  "The  Geometry  for  Oxetane  -  Protonated  Oxetane  Complex" 

The  optimal  angle  a  between  the  rings  has  been  found  to  be  90°,  but 
differences  between  different  angles  are  not  significant. 

The  opening  (C4A-01A)  of  the  protonated  oxetane  ring  starts  at 

R(0ib-C4A)  =*  4.6  bohrs.  Next  both:  0^ b— ^4A  anc^  ^4A~^1A  chan9e 

simultaneously  until  the  complex  reaches  the  stabilization  point  at 
R(0iB-C4A)  at  2.9  bohrs  and  6  -  19°  (fully  open).  The  stabilization  energy 

at  R  *  2.9  bohrs  and  6  -  19°  is  -.04378  a.u.  -  -27.97  kcal/mol.  The 
estimated  activation  energy  is  6.27  kcal/mol. 

b'.  Method  of  calculation 

The  ab-initio  M0DP0T/VRDD0 

calculations  have  been  carried  out  for  each  point  of  the  surface.  Because 
of  the  size  of  the  intermolecular  complex  of  the  two  reacting  species 
molecular  orbitals  selected  from  the  localized  space  were  used  in  the  MRD-CI 
calculations.  The  geometries  studied  include  the  most  sensitive  part  of 
the  complex  (Figure  1 1 -8 )  in  the  Cl  procedure.  Ten  of  the  most  important 
main  reference  configurations  have  been  used  in  MRD-CI  treatment,  and  the 
same  set  of  main  reference  configurations  has  been  kept  through  whole 
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potential  surface.  The  use  of  mul ticonf igurational  scheme  is  to  assure 
avoiding  of  the  possible  artifacts. 


The  potential  energy  surfaces  and  reaction  potential  map  of  OXET  - 


OXETH  complex  are  presented  in  Figures  1 1 -2  to  1 1 -5 . 


Figure  1 1 —2 :  MRD-CI  Extrapolated  Energy  for  Oxetane-Protonated  Complex  for 
Fixed  Angle  6  and  Different  Intermolecular  Distances  R ( 01 B-C4A) 

Figure  1 1 -3 :  MRD-CI  Extrapolated  Energy  for  Oxetane-Protonated  Oxetane 

Complex  for  Fixed  Intermolecular  Distances  R(01B-C4A)  and  Different  6 
Angle  Values 


Figure  1 1 -4 :  0XET-0XETH  ,  Extrapolated  Cl  Energy  Along  the  Reaction 
Coordinate  for  Oxetane  plus  Protonated  Oxetane  Addition  Reaction. 
Figure  1 1 -5 :  The  Potential  Energy  Surface  for  OXET  Approaching  Protonated 
OXET 


The  detailed  tables  of  results  follow. 


Table  1 1-1 :  "OXET  +  OXETH  [R(01 B-C4A)  -  3.4  bohrs]  (different  a 
values) 

Table  1 1 -2 :  OXET  +  0XETH+6  =  0°  (fully  closed).  Energies  (a.u.)  as  a 
function  of  R ( 01 B-C4 A) 

Table  1 1  -3 :  OXET  +  0XETH+<5  *  5°  ,  Energies  (a.u.)  as  a  function  of 
R(01B-C4A) 


Table  1 1 -4 :  OXET  +  OXETH  6-10°  ,  Energies  (a.u.)  as  a  function  of 
R  ( 01 B— C4A ) 

Table  1 1 -5 :  OXET  +  0XETH+6  -  15°  ,  Energies  (a.u.)  as  a  function  of 
R ( 01 B-C4A) 


Table  1 1 -6 :  OXET  +  OXETH  5  -  19°  ,  Energies  (a.u.)  as  a  function  of 
R  ( 01 B— C4  A ) 


C4A)  (bohr) 


REACTION  COORDINATE 


R(0IB-C4A )  (bohr) 


19 


TABLE  II  -  1 


OXET  +  OXETH*  [R(01 B-C4A)  -  3.4  bohrs] 
(DIFFERENT  ct  VALUES) 

ENERGIES  (a.u.) 


6=0°  (FULLY  CLOSED) 

6  -  19° 

(FULLY  OPEN) 

SCF 

Cl, EX 

SCF 

Cl, EX 

90 

135 

180 

-71.165534 

-71.163113 

-71.160545 

-71 . 359265 
-71.355541 
-71.352441 

-71.120867 

-71.120512 

-71.119567 

-71.306027 

-71.304423 

-71.301750 

minimum  for  a  -  90° 

a  =  angle  between  planes  of  rings 

<5  -  angle  of  opening  of  protonated  ring 


TABLE  II  -  2 


A 


OXET  +  OXETH+ 
<5=0°  (FULLY  CLOSED) 

ENERGIES  (a.u. ) 


R ( 01 8-C4 A) 


(bohrs) 

2.1 

2.5 

2.9 

SCF 

-70.560381 

-70.935433 

-71.133963 

Cl 

-70.755971 

-71.130522 

-71.325511 

EX 

-70.764130 

-71.138320 

-71.332365 

DAV 

-70.767576 

-71.141576 

-71.335368 

c2 

.963 

.964 

.964 

gs 

.904 

.905 

.905 

R( 01 B-C4A) 

(bohrs) 

3.6 

4.1 

4.6 

SCF 

-71.282489 

-71.3153506 

-71.3248406 

Cl 

-71.467731 

-71.498166 

-71.506188 

EX 

-71.474386 

-71.503304 

-71.511483 

DAV 

-71.476889 

-71.505569 

-71.513815 

c2 

.967 

.968 

.989 

gs 

.905 

.905 

.905 

R( 01 B-C4A) 

(bohrs) 

10.0 

SCF 

-71.315625 

Cl 

-71.497548 

EX 

-71.501413 

DAV 

-71.503505 

c2 

.969 

gs 

.903 

2 

c  is  the  contribution  of  all  of  the  reference  configurations 
gs  is  the  contribution  of  the  ground  state  SCF  wave  functions 


TABLE  II  -  3 


B 


OXET  +  OXETH+ 


(6 


) 


ENERGIES  (a.u.) 


R(01B-C4A) 
(bohrs ) 

2.1 

2.5 

2.9 

SCF 

-70.304448 

-71.109877 

-71.241954 

Cl 

-70.991765 

-71.298261 

-71.429245 

EX 

-70.998971 

-71.305820 

-71.437324 

DAV 

-71.001469 

-71.308416 

-71.440114 

? 

cu 

.968 

.967 

.965 

gs 

.905 

.905 

.904 

P( 01 3-C4A) 

(bohrs) 

3.6 

4.1 

4.6 

SCF 

-71.315482 

-71.324034 

-71.322182 

Cl 

-71.498182 

-71.505390 

-71.502834 

EX 

-71.503904 

-71.510932 

-71.507776 

DAV 

-71.506605 

-71.513540 

-71.510269 

2 

c 

.964 

.964 

.965 

gs 

.904 

.902 

.901 

R( 01 B-C4A) 

(bohrs ) 

10.0 

SCF 

-71.303100 

Cl 

-71.485062 

EX 

-71.488757 

DAV 

-71.491173 

c2 

.965 

gs 

.898 
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TABLE  II  -  4 
C 


R ( 01 B-C4A) 


OXET  +  OXETH 
(<5  =  10°) 

ENERGIES  (a.u.) 


(bohrs) 

2.1 

2.5 

2.9 

SCF 

-70.9586^ 

-71.217923 

-71.308731 

Cl 

-71.138404 

-71.399798 

-71.491715 

EX 

-71.144042 

-71.405662 

-71.497156 

DAV 

-71.146096 

-71.407782 

-71.499457 

c2 

.971 

.970 

.968 

gs 

.908 

.907 

.905 

R ( 0 1 B-C4A) 
(bohrs) 

3.6 

4.1 

4.6 

SCF 

-71.330412 

-71.318510 

-71.305275 

Cl 

-71.509340 

-71.494130 

-71.478736 

EX 

-71.514915 

-71.499211 

-71.483719 

DAV 

-71.517429 

-71.501685 

-71.486078 

2 

c 

.965 

.964 

.965 

gs 

.903 

.904 

.903 

R(0I B-C4A) 
(bohrs ) 

10.0 

SCF 

-71.273589 

Cl 

-71.446907 

EX 

-71.450690 

DAV 

-71.452995 

2 

c 

.964 

gs 

.901 

S’> 


»:•  i 


TABLE  II  -  5 


R(01B-C4A) 


R(01 B-C4A) 


R( 01 B-C4A) 
bohrs) 


-71.246902 

-71.414410 

-71.417747 

-71.419465 


OXET  +  OXETH 


(6  =  15u) 


ENERGIES  (a.u.) 


-71.052683 

-71.227256 

-71.232863 

-71.234740 


-71.338394 

-71.516580 

-71.521608 

-71.523932 


-71.281971 

-71.459679 

-71.464859 

-71.466807 


-71.347995 

-71.527670 

-71.532460 

-71.535503 


-71.312645 

-71.486829 

-71.491883 

-71.494220 


-71.290592 

-71.460933 

-71.465153 

-71.467242 
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TABLE  II  -  6 


E 

OXET  +  OXETH+ 

6  -  19°  (FULLY  OPEN) 

ENERGIES  (a.u.) 


R ( 01 B-C4A) 
(bohrs) 

2.1 

SCF 

-71.101569 

Cl 

-71.274272 

EX 

-71.279263 

DAV 

-71.281108 

c2 

.973 

gs 

.913 

R( 01 B-C4A) 

(bohrs) 

3.6 

SCF 

-71.332170 

Cl 

-71.512029 

EX 

-71.516097 

DAV 

-71.518387 

c2 

.966 

gs 

.901 

R( 01 B-C4A) 

(bohrs ) 

10.0 

SCF 

-71.217258 

Cl 

-71.383225 

EX 

-71.385572 

DAV 

-71.387004 

c2 

.973 

gs 

.905 

2.5 

2.9 

-71.310287 

-71.360621 

-71.486679 

-71.540063 

-71.491058 

-71.545193 

-71.492874 

-71.547200 

.972 

.970 

.910 

.906 

4.1 

4.6 

-71.297912 

-71.269955 

-71.474238 

-71.441685 

-71.478066 

-71.445050 

-71.480429 

-71.447196 

.964 

.966 

.900 

.902 

Ev 
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Oxetane  (OXET)  +  Protonated  3-fluoro-3- 
nitrooxetane  (FN0XH+) 

a*.  Determination  of  the  pathway  for 

oxetane  (OXET)  and  protonated  3-fluoro-3-ni trooxetane  (FN0XH+)  addition 
reaction . 


The  addition  of  oxetane  (OXET)  to 
protonated  3-fl uoro-3-ni trooxetane  (FN0XH+)  has  the  features  of  the  S^2 
reaction,  and  the  approach  presented  utilizes  these  to  the  greatest  extent. 

The  geometry  of  reacting  system. 

The  reaction  of  oxetane  with 

protonated  FNOX  is  assumed  to  be  the  S^2  type  reaction,  with  molecule  B 

attacking  molecule  A  (Figure  II  -  6)  along  the  (Lg-C^-O^  line. 

Figure  1 1 -6 :  "The  Geometry  For  Oxetane  -  Protonated  FNOX  Complex"  Since  it 
was  found  for  oxetane  reacting  with  protonated  oxetane  that  head  on  geometry 
for  approaching  oxetane  is  energetical ly  preferred  the  atom  is  also  kept 

on  °ib~('4A_01A  line.  Because  essential  structural  changes  are  expected  to 

appear  relative  to  the  plane  of  FNOX  ring  the  angle  alpha  between  the  two 
rings  was  determined  for  R(0  1B-C4A)  "  2’6  bohrs  and  was  used  for  all  other 

geometries.  The  optimal  alpha  angle  has  been  found  to  be  55  degrees. 

Figure  1 1 -7 ,  Table  1 1 —7  "OXET  +  FN0XH+,  SCF  And  Cl  Energy  For  Different 
Values  of  Angle  Between  Two  Rings  (Alpha)  (Fully  Closed  Geometry)" 

Two  bonds:  C4A-01A  and  01B-C4A  are  essential  to  describe  the  reaction 
pathway.  The  bond  inside  the  FNOX  ring  (C4A-01A)  changes  only  from  R  -  2.8 
bohrs  to  R  »  4.9  bohrs,  which  corresponds  to  the  fully  closed  and  fully  open 

ring  of  FN0XH+.  This  bond  is  later  described  by  the  more  natural  parameter 
6  which  varies  from  0°  (fully  closed)  to  19°  (fully  open),  and  corresponds 
to  degree  of  openess  of  the  ring.  The  01B-C4A  bond  is  changed  from  R  -  2.1 
bohrs  to  R  -  10.0  bohrs. 


Positions  of  proton  H+  and  hydrogens  connected  to  C4A  atom  are  the  most 
affected  by  opening  the  ring  and  these  positions  have  been  determined 
optimally  for  each  different  FN0XH+-ring  geometry. 

b'.  Method  of  Calculation 

The  ab-initio  M0DP0T/VRDD0 

calculations  have  been  carried  out  for  each  point  of  the  surface.  Because 
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of  the  size  of  the  i ntermol ecular  complex  of  the  two  reacting  species 
molecular  orbitals  selected  from  the  localized  space  were  used  in  the  MRD-CI 
calculations.  The  geometries  studied  include  the  most  sensitive  part  of  the 

complex  (Figure  II  -  8:  "The  Geometrical  State  of  Oxetane-FNOXH+  Selected 
for  MRO-CI  Calculations")  in  the  Cl  procedure.  Ten  of  the  most  important 
main  reference  configurations  have  been  used  in  MRD-CI  treatment,  and  the 
same  set  of  main  reference  configurations  has  been  kept  through  whole 
potential  surface.  The  use  of  mul ticonf igurational  scheme  is  to  assure 
avoiding  of  possible  artifacts. 


c1.  Results  for  oxetane  and  FN0XH+. 

The  ring  opening  starts  when 

oxetane  approaches  FN0XH+  for  R(Ojg-C^)  =  4.6  bohrs.  Next,  both: 

and  C4^-0jA  change  simultaneously  until  the  complex  reaches  the 

stabilization  point  at  -  2.9  bohrs.  and  6  -  19°  (fully  open). 

The  stabilization  energy  for  the  complex  is  -0.0603  a.u.  -  -39.7  kcal/mol. 

The  reaction  goes  through  transition  state  with  activation  energy 

approximately  .005  a.u.  =  3.13  kcal/mol. 

The  potential  energy  surfaces  and  reaction  potential  map  for  oxetane- 
+ 

FN0XH  complex  are  presented  in  Figures  1 1-9  to  1 1 -1 1 . 

Figure  1 1-9:  "MRD-CI  Extrapolated  Energy  for  Oxetane  Approaching  Protonated 
FN0X  for  Fixed  Angle  6  and  Different  Intermolecular  Distances 
R(01B-C4A) 

Figure  1 1 -1 0 :  "MRD-CI  Extrapolated  Energy  for  Oxetane-Protonated  FN0X 

Complex  For  Fixed  Intermolecular  Distances  R(01B-C4A)  and  Different  6 
Angle  Values 

Figure  1 1 - 1 1 :  "0XET  +  FN0XH  +  ,  Extrapolated  Cl  Energy  Along  The  Reaction 
Coordinate  For  Oxetane-Protonated  FN0X  Addition  Reaction 

Figure  11-12:  "The  Potential  Energy  for  Oxetane  Approaching  FN0X" 

The  detailed  tables  of  results  follow. 

Table  1 1-8 :  "0XET  +  FN0XH+  <5=0°  (fully  closed),  Energies  (a.u.)  as  a 
function  of  R(01B-C4A)" 

Table  1 1 —9 :  "0XET  +  FN0XH+  6  =  5°,  Energies  (a.u.)  as  a  function  of 
R( 01B-C4A) " 

Table  1 1  - 1 0 :  "OXET  +  FN0XH+  <5  =  10°,  Energies  (a.u.)  as  a  function  of 
R ( 01 B-C4A) " 

Table  1 1 - 1 1 :  "OXET  +  FN0XH+  <5  =  15°,  Energies  (a.u.)  as  a  function  of 
R ( 01 B-C4A) " 


Table  11-12:  "OXET  +  FNOXH  <5  =  19  (fully  open),  Energies  (a.u.)  as  a 
function  of  R(01B-C4A)" 
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R(0IB-C4A)  (bohr) 


REACTION  COORDINATE 


THE  POTENTIAL  ENERGY  SURFACE  FOR  OXETANE  APPROACHING  PROTONATED  FNOX. 


RI0IB-C4A)  (bohr) 


TABLE  II  -  7 


OXET  +  FNQXH+ 

SCF  AND  Cl  ENERGY  (a.u.)  FOR  DIFFERENT  VALUES 
OF  ANGLE  BETWEEN  TWO  RINGS  (ALPHA) 

(FULLY  CLOSED  GEOMETRY) 


ot  (°) 

-0 

-45 

-55 

SCF 

-134.40904 

-134.56213 

-134.56319 

Cl 

-134.59884 

-134.75371 

-134.75493 

EX 

-134.60698 

-134.76151 

-134.76290 

DAV 

-134.61009 

-134.76467 

-134.76608 

,o. 

a  (  ) 

-65 

90  (-90) 

-135 

SCF 

-134.55734 

-134.43763 

-134.25124 

Cl 

-134.74906 

-134.62844 

-134.43975 

EX 

-134.75688 

-134.63601 

-134.44762 

DAV 

-134.76006 

-134.63916 

-134.45070 

TABLE  II  - 

A 

R ( 01 B-C4A) 
(bohrs ) 

8 

OXET  +  FN0XH+ 

6  =  0°  (FULLY  CLOSED) 

ENERGIES  (a . u . ) 

2.1  2.6 

3.1 

SCF 

-134.110905 

-134.563192 

-134.762800 

Cl 

-134.302148 

-134.754932 

-134.950999 

EX 

-  1 34. 310289 

-134.762896 

-134.958205 

DAV 

-134.313487 

-134.766081 

-134.961110 

2 

c 

.965 

.964 

.965 

GS 

.907 

.907 

.907 

R  ( 01 B-C4A) 

(bohrs ) 

3.6 

4.1 

4.6 

SCF 

-134.847935 

-134.878792 

-134.887225 

Cl 

-135.031883 

-135.060623 

-135.067424 

EX 

-135.038693 

-135.065956 

-135.072717 

OAV 

-135.041224 

-135.068729 

-135.074857 

2 

c 

.967 

.962 

.969 

gs 

.907 

.902 

.906 

R(01 B-C4A) 

(bohrs ) 

7.1 

10.0 

SCF 

-134.882235 

-134.877230 

Cl 

-135.062520 

-135.058033 

EX 

-135.066463 

-135.061843 

DAV 

-135.068503 

-135.063898 

2 

c 

.970 

.970 

gs 

.905 

.904 

vLvALvy.V'. 
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TABLE  II  -  9 
B 

OXET  +  FNOXH+ 
6  =  5° 


ENERGIES  (a.u.) 


R(01B-C4A) 
(bohrs ) 

2. 1 

2.6 

3.1 

SCF 

-134.370920 

-134.727815 

-134.847435 

Cl 

-134.553873 

-134.914168 

-135.033105 

EX 

-134.560230 

-134.921401 

-135.040733 

DAV 

-134.563793 

-134.932963 

-135.043589 

c2 

.961 

.967 

.964 

gs 

.908 

.906 

.905 

R( 01 B-C4A) 

(bohrs ) 

3.6 

4.1 

4.6 

SCF 

-134.882307 

-134.887248 

-134.883394 

Cl 

-135.065067 

-135.068662 

-135.063898 

EX 

-135.070494 

-135.074395 

-135.068526 

DAV 

-135.074265 

-135.077118 

-135.071054 

c2 

.957 

.964 

.965 

gs 

.904 

.903 

.901 

R( 01 B-C4A) 

(bohrs ) 

7. 1 

10.0 

SCF 

-134.868181 

-134.862798 

Cl 

-135.049391 

-135.044737 

EX 

-135.052791 

-135.048036 

DAV 

-135.055166 

-135.050435 

c2 

.965 

.965 

gs 

.899 

.898 

TABLE  II  -  10 


C 

OXET  ♦  FNOXH+ 
6  -  10° 

ENERGIES  (a.u.) 


R(01B-C4A) 
(bohrs ) 

2.1 

2.6 

3.1 

SCF 

-134.535644 

-134.830287 

-134.899875 

Cl 

-134.711865 

-135.010989 

-135.082938 

EX 

-134.717174 

-135.016806 

-135.088818 

DAV 

-134.719132 

-135.018927 

-135.091260 

2 

c 

.971 

.970 

.967 

gs 

.910 

.907 

.904 

R( 01 B-C4A) 

(bohrs ) 

3.6 

4.1 

4.6 

SCF 

-134.899562 

-134.882823 

-134.866316 

Cl 

-135.080986 

-135.060244 

-135.040898 

EX 

-135.086011 

-135.065866 

-135.045438 

DAV 

-135.088703 

-135.068135 

-135.048004 

2 

c 

.964 

.968 

.963 

gs 

.901 

.906 

.903 

R( 01 B-C4A) 

(bohrs) 

7.1 

10.0 

SCF 

-134.837119 

-134.831005 

Cl 

-135.010516 

-135.005397 

EX 

-135.014918 

-135.009072 

DAV 

-135.017270 

-135.011458 

2 

c 

.964 

.967 

g? 

.900 

.900 

TABLE  II  -  11 


D 

OXET  +  FNOXH+ 

6  -  15° 

ENERGIES  (a.u.) 


R ( G1 B-C4 A) 
(bohrs ) 

2.1 

2.6 

3.1 

SCF 

-134.635098 

-134.889992 

-134.930457 

Cl 

-134.806808 

-135.066946 

-135.111510 

EX 

-134.811851 

-135.072539 

-135.117082 

DAV 

-134.813663 

-135.074458 

-135.119276 

c2 

.973 

.972 

.968 

gs 

.913 

.909 

.904 

R(01B-C4A) 

(bohrs) 

3.6 

4.1 

4.6 

SCF 

-134.909566 

-134.878881 

-134.853019 

Cl 

-135.091031 

-135.057257 

-134.026524 

EX 

-135.096412 

-135.062585 

-134.030965 

DAV 

-135.098936 

-135.065315 

-134.033521 

2 

c 

.965 

.962 

.962 

GS 

.899 

.898 

.900 

R( 01 B-C4A) 
(bohrs) 

7.1 

7.1 

SCF 

-134.810279 

-134.810279 

Cl 

-135.977164 

-134.977164 

EX 

-135.980338 

-134.980388 

DAV 

-135.982137 

-134.982137 

2 

c 

.969 

.969 

gs 

.904 

.904 

m 


i 


W' 
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TABLE  I!  -  12 


OXET  +  FNOXH 
5  =  19°( FULLY  OPEN) 

ENERGIES  (a.u.) 


R( 01 B-C4A) 
(bohrs ) 

2.1 

2.6 

2.9 

SCF 

-134.685867 

-134.914697 

-134.940369 

Cl 

-134.856194 

-135.090774 

-135.119754 

EX 

-134.861224 

-135.095386 

-135. 125109 

DAV 

-134.863028 

-135.097236 

-135.127119 

c2 

.973 

.972 

.970 

gs 

.915 

.910 

.907 

R(01B-C4A) 

(bohrs) 

3.1 

3.6 

4.6 

SCF 

-134.937775 

-134.905172 

-134.866571 

Cl 

-135.118956 

-135.088558 

-135.048086 

EX 

-135. 123972 

-135.092678 

-135.052252 

DAV 

-135.126087 

-135.095123 

-135.055025 

c2 

.969 

.965 

.961 

gs 

.904 

.897 

.894 

R( 01 B-C4A) 

(bohrs ) 

4.6 

7.1 

10.0 

SCF 

-134.835066 

-134.782966 

-134.775784 

Cl 

-135.011475 

-134.948137 

-135.942239 

EX 

-135.014585 

-134.951378 

-135.945234 

DAV 

-135.017295 

-134.952855 

-135.946753 

c2 

.960 

.973 

.972 

gs 

.895 

.905 

.904 
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TABLE  11-13 


a 


3 


m 


t 


FNOX  +  OXETH 


6  =  Ou  (FULLY  CLOSED) 


SCF  ENERGY  (a.u.)  FOR  DIFFERENT  VALUES  OF  ANGLE  BETWEEN  TWO  RINGS  (a) 


«  n 


-134.573887 

-134.566354 

-134.556369 

-134.566333 

-134.573804 
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3‘.  3-F 1 uoro-3-ni trooxetane  (FNOX)  +  protonated 

oxetane  (0XETH+). 

a'.  Determination  of  the  pathway  for  3- 
fluoro-3-nitrooxetane  (FNOX)  and 

protonated  oxetane  (0XETH+)  addition. 
The  addition  of  3-fluoro-3- 

nitrooxetane  (FNOX)  to  protonated  oxetane  (0XETH+)  also  has  the  features  of 
an  S^2  reaction  and  the  approach  prescibed  utilizes  these  to  the  greatest 

extent. 

The  geometry  of  the  reacting  system. 

The  reaction  of  FNOX  with  OXETH*  is  assumed  to  be  SN2  type  reaction 

with  molecule  B  attacking  molecule  A  (Figure  11-13)  along  the  °ib_^4A_^1A 

line.  Because  essential  structural  changes  are  expected  to  appear  relative 
to  the  plane  of  FNOX  ring  the  angle  alpha  between  the  two  rings  was 
determined  for  R(01B-C4A)  -  2.6  bohrs  and  was  used  for  all  other  geometries. 

b'.  Method  of  calculation 

There  is  little  difference  in 

energy  with  changes  in  a,  the  angle  between  the  rings,  Table  11-13,  "FNOX  + 

0XETH+,  6=0°  (fully  closed),  SCF  Energy  for  Different  Values  of  Angle 
Between  Two  Rings  (a)".  The  method  of  calculation  is  the  same  as  described 

for  OXET  +  FN0XH+  in  the  previous  section. 

c'.  Results  for  FNOX  +  0XETH+. 

The  ring  opening  starts  when  FNOX 

approaches  to  oxetane  H+  for  R(0^g-C^^)  =  4.6  bohrs.  Next,  both:  ^b-^ 
and  C4A~0 ^  change  simultaneously  until  the  complex  reaches  the 

stabilization  point  at  *  2.9  bohrs  and  6  -  19°  (fully  open).  The 

stabi 1  ization  energy  for  the  complex  is  -.00915  a.u.  -  -5.74  kcal/mol.  The 
estimated  activation  energy  is  6.27  kcal/mol. 

The  energies  for  angle  a  between  rings  do  not  vary  significantly. 

The  potential  energy  surfaces  and  reaction  potential  map  are  presented 
in  Figures  11-14  to  Figure  11-17. 


Figure  11-14:  "MRD-CI  Extrapolated  Energy  For  FNOX  Approaching  Protonated 
Oxetane  For  Fixed  Angle  6  And  Different  I ntermolecul ar  Distances 
R ( 01 S-C4A) " 

Figure  11-15:  "MRD-CI  Extrapolated  Energy  For  FNOX-Protonated  Oxetane 

Complex  For  Fixed  Intermolecular  Distances  R(01B-C4A)  and  Different  6 
Angle  Values" 

Figure  11-16:  "FN0X-0XETH+  Extrapolated  Cl  Energy  Along  the  Reaction 
Coordinate  For  FNOX-protonated  Oxetane  Addition  Reaction" 

Figure  11-17:  "The  Potential  Energy  Surface  For  FNOX  Approaching  Protonated 
Oxetane" 

The  detailed  Tables  of  results  follow: 

Table  11-14:  "FNOX  +  0XETH+,  6=0°  (fully  closed)  Energies  (a.u.)  as  a 
function  of  R(01B-C4A)" 

Table  11-15:  "FNOX  +  OXETH*,  6=5°  Energies  (a.u.)  as  a 
function  of  R(01B-C4A)" 

Table  11-16:  "FNOX  +  0XETH+,  6  =  10°  Energies  (a.u.)  as  a 
function  of  R(01B-C4A)" 

Table  11-17:  "FNOX  +  0XETH+,  6  =  15°  Energies  (a.u.)  as  a 
function  of  R(01B-C4A)" 

Table  11-18:  "FNOX  +  0XETH+,  6  -  19°  (fully  open)  Energies  (a.u.)  as  a 
function  of  R(01B-C4A)" 
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TABLE  II  -  13 


FNOX  +  OXETH+ 
6=0°  (FULLY  CLOSED) 


SCF  ENERGY  (a.u.)  FOR  DIFFERENT  VALUES  OF  ANGLE  BETWEEN  TWO  RINGS  (a) 


SCF 

-134.573887 

-134.566354 

-134.556369 

-134.566333 

-134.573804 


TABLE  11-14 
A 

FNOX  +  OXETH+ 

6  =  0°  (FULLY  CLOSED) 

ENERGIES  (a.u.) 


R(01B-C4A) 

(bohrs) 

2.1 

2.5 

2.9 

SCF 

-134.131955 

-134.511150 

-134.71439 

Cl 

-134.328233 

-134.706202 

-134.905989 

EX 

-134.336047 

-134.713390 

-134.912682 

DAV 

-134.339619 

-134.716688 

-134.915661 

2 

c 

.962 

.963 

.964 

gs 

.904 

.905 

.905 

R(01B-C4A) 

(bohrs) 

3.6 

4.1 

4.6 

SCF 

-134.870168 

-134.906321 

-134.918159 

Cl 

-135.055303 

-135.088869 

-135.099405 

EX 

-135.061614 

-135.094212 

-135.104669 

DAV 

-135.064059 

-135.096435 

-135.106795 

2 

c 

.967 

.968 

.969 

gs 

.905 

.905 

.905 

R( 01 B-C4A) 

(bohrs) 

10.0 

SCF 

-134.917271 

Cl 

-135.099207 

EX 

-135.103049 

DAV 

-135.105145 

c2 

.969 

gs 

.903 
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TABLE  11-15 


FNOX  +  0XETH+ 
<5  =  5° 

ENERGIES  (a.u.) 

R ( 01 B-C4A) 


(bohrs ) 

2. 1 

2.9 

3.6 

SCF 

-134.371600 

-134.817032 

-134.899072 

Cl 

-134.559469 

-135.003579 

-135.081621 

EX 

-134.566346 

-135.010938 

-135.087747 

DAV 

-134.568878 

-135.013692 

-135.090412 

c2 

.968 

.965 

.964 

gs 

.905 

.905 

.904 

R (01 B-C4A) 

(bohrs ) 

4.1 

10.0 

SCF 

-134.912246 

-134.904573 

Cl 

-135.093753 

-135.086645 

EX 

-135.098762 

-135.090328 

DAV 

-135.101315 

-135.092754 

2 

c 

.965 

.965 

gs 

.902 

.898 
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TABLE  11-16 


R( 01 B-C4A) 
bohrs ) 


FNOX  +  OXETH 


6  =  10L 


ENERGIES  (a.u.) 


R(01B-C4A) 
(bohrs ) 

2.1 

2.5 

2.9 

SCF 

-134.523543 

-134.785323 

-134.879961 

Cl 

-134.703562 

-134.967402 

-135.062094 

EX 

-134.709561 

-134.973213 

-135.067821 

OAV 

-134.711637 

-134.875357 

-135.070119 

c2 

.971 

.970 

.968 

gs 

.908 

.907 

.905 

R( 01 3-C4A) 

(bohrs ) 

3.6 

4.1 

4.6 

SCF 

-134.909633 

-134.903107 

-134.894288 

Cl 

-135.087083 

-135.077655 

-135.067303 

EX 

-135.092434 

-135.082991 

-135.072249 

DAV 

-135.094867 

-135.085375 

-135.074560 

2 

c 

.965 

.965 

.964 

gs 

.905 

.904 

.901 

-134.874769 

-135.048181 

-135.051947 

-135.054272 
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TABLE  11-17 


FNOX  +  OXETH 
6  =  15° 

ENERGIES  (a.u.) 


R(01B-C4A) 


-134.616620 

-134.791828 

-134.797443 

-134.799299 


-134.917116 

-135.096343 

-135.101219 

-135.103258 


-134.914524 

-135.090622 

-135.096394 

-135.098661 


R( 01 B-C4A) 


-134.89416562 

-135.066025 

-135.070201 

-135.072315 


-134.847906 

-135.015423 

-135.018739 

-135.020467 


TABLE  11-18 


FNOX  +  OXETH 
6  =  19°  (FULLY  OPEN) 

ENERGIES  (a.u.) 


R ( 01 B-C4A) 
(bohrs ) 

2.1 

2.5 

2.9 

SCF 

-134.665293 

-134.875744 

-134.928616 

Cl 

-134.838385 

-135.052121 

-135.107297 

EX 

-134.843237 

-135.056507 

-135.112204 

DAV 

-134.845089 

-135.058348 

-135.114180 

c2 

.973 

.972 

.970 

gs 

.913 

.910 

.907 

R(01 B-C4A) 

(bohrs ) 

3.6 

4.6 

10.0 

SCF 

-134.906366 

-134.854626 

-134.818210 

Cl 

-135.083692 

-135.023718 

-134.984149 

EX 

-135.087620 

-135.027136 

-134.986463 

OAV 

-135.089794 

-135.028994 

-134.987900 

c2 

.967 

.969 

.973 

gs 

.903 

.904 

.904 
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4'.  3-Fluoro-3-ni trooxetane  (FNOX)  and 
protonated  3-f 1 uoro-3-ni trooxetane 

( FN0XH+) 

a'.  Determination  of  the  pathway  for 
3-fluoro-3-nitrooxetane  (FN0XH+) 

The  addition  of  3-fluoro-3- 

nitrooxetane  (FNOX)  to  protonated  3-fluoro-3-nitrooxetane  (FN0XH+)  also  has 
the  features  of  an  S^2  reaction  and  the  approach  presented  utilizes  these  to 

the  greatest  extent. 


The  geometry  of  the  reacting  system. 

The  reaction  of  FNOX  with  FN0XH+  is  assumed  to  be  S^2  type  reaction  with 

molecule  B  attacking  molecule  A  (Figure  11-18)  along  the  O^b'C^-O^  line. 

Because  essential  structural  changes  are  expected  to  appear  relative  to  the 
plane  of  FNOX  ring  the  angle  alpha  between  the  two  rings  was  determined  for 
R(01B-C4A)  =  2.6  bohrs  and  was  used  for  all  other  geometries. 

The  optimal  a  angle  was  determined  to  be  55°  (Table  11-19),  "FNOX  + 
FN0XH+,  SCF  Energy  for  Various  Values  of  a  Angle". 

b'.  Method  of  calculation 

The  method  of  calculation  is  the 
same  as  described  for  OXET  +  FN0XH+  in  section  V. 

c'.  Results  for  FNOX  +  FN0XH+ 

The  ring  opening  again  starts 

when  FNOX  approaches  FNOXH*  for  R( 01 B-C4 A)  =  4.6  bohrs.  Again  next  both: 

0 1 B_C 4 a  ancl  ^4A_01A  chan9e  simultaneously  until  the  stabilization  point  at 

R ^ 0 1 b_C4 A ^  =  2,9  bohrs  ar,d  <5  =  I90  (fully  open) 


The  stabi 1  ization  energy  for  the  complex  is  -.03156  a.u.  =  -19.80 
kcal/mol.  Approximate  activation  energy  .005  a.u.  =■  3.13  kcal/mol.  The 
potential  energy  surfaces  and  reaction  potential  map  are  presented  in 
Figures  11-19  to  11-22. 

Figure  11-19:  "MRD-CI  Extrapolated  Energy  For  FNOX  Approaching  Protonated 
FNOX  For  Fixed  Angle  6  and  Different  Intermolecular  Distances  R( 01 B— 
C4A) " 

Figure  11-20:  "MRD-CI  Extrapolated  Energy  For  FN0X-FN0X  Protonated  Complex 
For  Fixed  Intermolecular  Distances  R( 01 B-C4 A)  and  Different  6  Angle 


Values" 


Figure  11-21:  "FNOX-FNOXH  ,  Extrapolated  Cl  Energy  Along  the  Reaction  Coor¬ 
dinate  for  FNOX  Protonated  FNOX  Addition  Reaction" 

Figure  11-22:  "The  Potential  Energy  Surface  for  FNOX  Approaching  Protonated 
FNOX" 

The  detailed  Tables  of  results  follow: 


Table  11-20:  "FNOX  +  FN0XH+,  6=0°  (fully  closed)  Energies  (a.u.)  as  a 
function  of  R ( 01 B-C4A) " 

Table  11-21:  "FNOX  +  FN0XH+,  6=5°  Energies  (a.u.)  as  a 
function  of  R(01B-C4A)" 

Table  11-22:  "FNOX  +  FN0X«\  6  =  10°  Energies  (a.u.)  as  a 
function  of  R(01B-C4A)" 

Table  11-23:  "FNOX  +  FN0XH+,  6  =  15°  Energies  (a.u.)  as  a 
function  of  R(01B-C4A)" 

Table  11-24:  "FNOX  +  FN0XH+,  6  =  19°  (fully  open)  Energies  (a.u.)  as  a 
function  of  R(01B-C4A)“ 
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TABLE  11-19 


FNOX  +  FNOXH+ 

SCF  ENERGY  (a.u.)  FOR  VARIOUS  VALUES  OF  a  ANGLE 


,0, 
a  {  ) 

ENERGY 

0 

-198.181177 

45 

-198.281031 

55 

-198.281216 

65 

-198.276447 

90 

-198.199361 

135 

-198.079309 

180 

-198.182049 

225 

-198.280452 

235 

-198.280818 

245 

-198.276348 

270 

-198.199587 

MINIMUM 
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TABLE  11-20 


FNOX  +  FNOXH 
6=0°  (FULLY  CLOSED) 

ENERGIES  (a.u.) 


R (01 B-C4A) 

(bohrs ) 

2.3 

2.9 

3.6 

SCF 

-197.910863 

-198.281216 

-198.435192 

Cl 

-198.103109 

-198.470845 

-198.618954 

EX 

-198.111535 

-198.477910 

-198.625229 

DAV 

-198.114770 

-198.480888 

-198.627675 

c2 

.964 

.964 

.967 

gs 

.906 

.906 

.906 

R  ( 0 1 B-C4  A) 

(bohrs) 

4.6 

7.1 

10.0 

SCF 

-198.480690 

-198.418577 

-198.479303 

Cl 

-198.634550 

-198.661953 

-198.656760 

EX 

-198.638251 

-198.665930 

-198.660445 

DAV 

-198.639642 

-198.667978 

-198.662301 

c2 

.974 

.969 

.971 

gs 

.915 

.904 

.906 
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TABLE  11-21 


FNOX  +  FNOXH 
6  =  5° 

ENERGIES  (a.u.) 


R(01 B-C4A) 
(bohrs ) 

2.3 

2.9 

3.6 

SCF 

-198. 128612 

-198.389364 

-198.465030 

Cl 

-198.313697 

-198.575069 

-198.647397 

EX 

-198.320611 

-198.582967 

-198.653065 

DAV 

-198.323023 

-198.585722 

-198.655796 

2 

c 

.968 

.965 

.964 

gs 

.907 

.906 

.904 

R ( 01 B-C4A) 

(bohrs ) 

4.6 

10.0 

SCF 

-198.474937 

-198.464728 

Cl 

-198.655547 

-198.646778 

EX 

-198.660043 

-198.650079 

DAV 

-198.662521 

-198.652485 

2 

c 

.965 

.965 

gs 

.901 

.898 
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TABLE  11-22 


C 

FNOX  +  FN0XH+ 
6  =  10° 

ENERGIES  (a.u.) 


R(01B-C4A) 


(bohrs ) 

2.3 

2.9 

3.6 

SCF 

-198.265018 

-198.456342 

-198.477340 

Cl 

-198.443306 

-198.638345 

-198.656431 

EX 

-198.448377 

-198.644268 

-198.661916 

DAV 

-198.450357 

-198.646572 

-198.664543 

c2 

.971 

.968 

.934 

gs 

.909 

.905 

.903 

R( 01 B-C4A) 

(bohrs ) 

4.6 

10.0 

SCF 

-198.454554 

-198.432698 

Cl 

-198.628227 

-198.607392 

EX 

-198.633078 

-198.611131 

DAV 

-198.635537 

-198.613549 

c2 

.964 

.963 

gs 

.903 

.899 

TABLE  11-23 


FNOX  +  FNOXH 
5  =  15° 

ENERGIES  (a.u.) 


R( 01 B-C4A) 
(bohrs ) 

2.3 

2.9 

3.6 

SCF 

-198.346074 

-198.495395 

-198.483930 

Cl 

-198.520243 

-198.674606 

-198.662553 

EX 

-198.525408 

-198.679879 

-198.667532 

DAV 

-198.527228 

-198.681927 

c2 

.972 

.970 

.966 

95 

.912 

.906 

.902 

R( C4A-01 B) 

(bohrs) 

4.6 

10.0 

SCF 

-198.437954 

-198.404983 

Cl 

-198.608248 

-198.569461 

EX 

-198.612647 

-198.572663 

DAV 

-198.614829 

-198.574290 

TABLE  II  -  24 


FNOX  +  FNOXH 
6  =  19°  (FULLY  OPEN) 

ENERGIES  (a.u.) 


R( 01 B-C4A) 
(bohrs ) 

2.3 

2.9 

3.6 

SCF 

-198.384740 

-198.507770 

-198.477497 

Cl 

-198.557601 

-198.686766 

-198.657925 

EX 

-198.562331 

-198.692012 

-198.661934 

DAV 

-198.564119 

-198.694011 

-198.664266 

2 

c 

.973 

.970 

.966 

gs 

.913 

.906 

.900 

R ( 0 1 B-C4A) 

(bohrs ) 

4.6 

10.0 

SCF 

-198.417703 

-198.377272 

Cl 

-198.589036 

-198.540339 

EX 

-198.592776 

-198.543087 

DAV 

-198.594981 

-198.544440 

c2 

.965 

.973 

gs 

.902 

.906 
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(b).  Population  Analyses 


Gerry  Manser  had  expressed  considerable 
interest  in  how  the  charges  (corresponding  to  the  gross  atomic  populations) 
on  the  0^  (oxygen  of  protonated  oxetane  ring),  (the  a  carbon  of  the 

protonated  oxetane  ring)  and  O^g  (oxygen  of  the  oxetane  ring)  varied  as  a 

function  of  substituent  and  reaction  pathway. 

In  last  year's  ONR  Annual  Report  1986,  Table  III  Page  24,  showed  that 
that  as  oxetane  and  protonated  oxetane  approached  each  other  that  the  intra¬ 
ring  TOP  of  the  C^-0jA  in  the  protonated  ring  got  smaller  as  the  oxetane 

ring  approached,  indicating  a  tendency  for  the  protonated  ring  to  open  and 
the  inter-ring  TOP  got  larger  indicating  bond  formation. 

On  the  following  pages  we  have  Tables  11-25  to  11-28  of  the  MRD-CI 
gross  atomic  populations  and  total  overlap  populations  (C^-O^.)  and 

(01B_C4A)  for  the  systems  0XET  +  0XETH+,  OXET  +  FN0XH+,  FNOX  +  0XETH+  and 

FNOX  +  FN0XH+  at  representat i ve  intermediate  inter-ring  distance  0^B~C4A 
distance  of  3.6  bohrs. 

i)  It  is  apparent  from  the  TOP'S  in  the  tables  that  the  two  rings  are 
repulsive  when  the  protonated  oxetane  (or  substituted  protonated  oxetane) 
ring  is  closed. 

ii)  the  protonated  oxetane  (or  substituted  protonated  oxetane)  will  open 
upon  approach  of  the  oxetane  (or  substituted  oxetane)  along  the  appropriate 
reaction  pathway. 

iff)  The  01B-C4A  interring  bond  becomes  stronger  as  the  protonated  (A) 
ring  opens. 

The  behavior  of  these  TOP'S  is  indicative  of  the  same  conclusion  as 
that  from  the  MRD-CI  energy  calculations. 
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(2).  AMMO  +  OXETH+,  OXET  ♦  AMMOH+,  AMMO  ♦  AMMOH+ 

The  reactions  AMMO  +  OXETH+,  OXET  +  AMMOH+  and 

AMMO  +  AMMOH+  were  treated  as  S^2  reactions  as  the  previous  cases  described 

in  detail  in  the  earl  ier  section.  The  optimal  a  angles  were  determined  to 
be  90°  for  each  of  those  systems  involving  AMMO  (3-azidomethyl-3- 

methy loxetane  and/or  protonated  AMMO  (AMM0Hf).  The  subsequent  MRD-CI 
calculations  were  run  at  this  a  angle. 

The  same  localized  skeletal  molecular  orbitals  (both  occupied  and 
virtual)  were  included  explicitly  in  the  MRD-CI  as  for  the  previous  cases. 
This  approach  makes  it  meaningful  to  compare  reaction  energies  between 
various  systems. 

For  these  cases  involving  AMMO  and/or  AMM0H+,  the  ring  opening  of  the 
protonated  ring  (0XETH+  or  AMM0H+)  starts  for  R(Ojg-C^)  -  4.6  bohrs.  Next, 
both:  0ig-C4A  and  C4A-01A  change  simultaneously  until  the  stabilization 
point  at  R(0^-C4A)  »  2.9  bohrs  and  6  =»  19°  (fully  open). 

(a).  Energies 

1'.  3-azidomethy 1-3-methyloxetane  (AMMO) 

+  protonated  oxetane  (0XETH+) 

a'.  Results 

The  stabilization  point  R( 01 B- 

C4A)  =  2.9  bohrs,  6  -  19°  (fully  open)  with  the  stabilization  energy 
-0.043262  a.u.  =  -27.37  kcal/mole. 

The  potential  energy  surfaces  and  reaction  potential  map  are  presented 
in  Figures  11-23  to  11-26. 

Figure  11-23:  "MRD-CI  Extrapolated  Energy  for  AMMO  Approaching  Protonated 
Oxetane  For  Fixed  Angle  6  and  Different  Intermolecular  Distances 
R ( 01 8-C4A) " 

Figure  11-24:  "MRD-CI  Extrapolated  Energy  For  AMMO-Protonated  Oxetane 

Complex  for  Fixed  Intermolecular  Distances  R ( 01 B-C4 A)  and  Different  5 
Angle  Values" 

Figure  11-25:  "AMMO  +  0XETH+,  Extrapolated  Cl  Energy  Along  the  Reaction 
Coordinate  For  AMMO  Protonated  Oxetane  Addition  Reaction." 

Figure  11-26:  "The  Potential  Energy  Surface  For  AMMO  Approaching  Protonated 
OXET,  MRD-CI  Extrapolated” 


Tne  detailed  Tables  of  results  follow  in  Tables  11-29  -  11-33: 


Table  11-29:  "AMMO  +  0XETH+  6  =  0°  (fully  closed),  Energies  (a.u.)  as  a 
function  of  R(01B-C4A)" 

Table  11-30:  "AMMO  +  OXETH*  6  -  5°,  Energies  (a.u.)  as  a  function  of 
R ( 01 B-C4A) " 

Table  11-31:  "AMMO  +  OXETH'*'  6  =  10°,  Energies  (a.u.)  as  a  function  of 
R ( 01 8-C4 A ) " 

Table  11-32:  "AMMO  +  0XETH+  6  =  15°,  Energies  (a.u.)  as  a  function  of 
R ( 01 8-C4A) " 

Table  11-33:  "AMMO  +  0XETH+  6  -  19°  (fully  open),  Energies  (a.u.)  as  a 
function  of  R(01B-C4A)" 


AMMO  ♦  OXETANEH  . EXTRAPOLATED  Cl  ENERGY  ALONG  THE  REACTION 


R(OlB-C4A)[bohr] 


R( 01B-C4A) 
(bohrs ) 

AMMO  +  OXETH 

6  -  0°  (FULLY  CLOSED) 

ENERGIES  (a.u.) 

2.3  2.9 

3.6 

SCF 

-112.212402 

-112.567419 

-112.715933 

Cl 

-112.407367 

-112.758700 

-112.901015 

EX 

-112.415741 

-112.765499 

-112.907699 

DAV 

-112.419053 

-112.768486 

-112.910196 

2 

c 

.964 

.964 

.967 

GS 

.905 

.906 

.905 

R(01 B-C4A) 

(bohrs) 

4.6 

10.0 

SCF 

-112.758393 

-112.750039 

Cl 

-112.939447 

-112.931844 

EX 

-112.944701 

-112.935694 

DAV 

-112.946828 

-112.937779 

2 

c 

.969 

.969 

gs 

.905 

.903 

rwrwi.rw' 


AMMO  + 

OXETH 

6  - 

5° 

ENERGIES 

(a.u.) 

R (01 S-C4A) 

(bohrs ) 

2.3 

2.9 

3 . 6 

SCF 

-112.422599 

-112.675770 

-112.748987 

Cl 

-112.610259 

-112.862895 

-112.931570 

EX 

-112.617265 

-112.871022 

-112.937359 

DAV 

-112.619748 

-112.873804 

-112.940058 

c2 

.968 

.966 

.965 

gs 

.906 

.905 

.905 

R( 01 B-C4A) 

10.0 

(bohrs) 

4.6 

SCF 

-112.755698 

-112.737486 

Cl 

-112.936221 

-112.919347 

EX 

-112.941127 

-112.923042 

DAV 

-112.943611 

-112.925452 

AMMO  + 

OXETH 

5  = 

10° 

ENERGIES 

;  (a.u.) 

R{ 01 B-C4A) 
(bohrs ) 

2.3 

2.9 

3.6 

SCF 

-112.553838 

-112.742890 

-112.764062 

Cl 

-112.734576 

-112.925800 

-112.943298 

EX 

-112.740024 

-112.931127 

-112.949404 

DAV 

-112.742076 

-112.933420 

-112.951970 

c2 

.971 

.968 

.965 

gs 

.908 

.905 

.904 

R ( 01 B-C4A) 
bohrs) 

4.6 

10.0 

— 1 - <— 

SCF 

-112.707936 

-112.738741 

Cl 

-112.912073 

-112.881148 

EX 

-112.917085 

-112.884915 

DAV 

-112.919438 

-112.887215 

AMMO  + 

OXETH 

6  = 

15° 

ENERGIES 

(a.u.) 

R ( 01 B-C4A) 
(bohrs ) 

2.3 

2.9 

3.6 

SCF 

-112.632473 

-112.782382 

-112.772218 

Cl 

-112.808287 

-112.962256 

-112.950352 

EX 

-112.813023 

-112.967384 

-112.955341 

DAV 

-112.814859 

-112.969441 

-112.957658 

2 

c 

.972 

.970 

.966 

gs 

.911 

.906 

.903 

R( 01 8-C4A) 
(bohrs) 

4.6 

10.0 

SCF 

-112.724061 

-112.632473 

Cl 

-112.894278 

-112.808287 

EX 

-112.898511 

-112.813023 

DAV 

-112.900595 

-112.814859 

80 


tabl:  : : -33 

r* 

AMMO  +  OXETH+ 

6  =  19°( FULLY  OPEN) 

ENERGIES  (a.u.) 


R ( 01 B-C4A) 

( bohrs ) 

2.3 

2.9 

3.6 

SCF 

Cl 

EX 

DAV 

-112.6704978 

-112.844764 

-112.849093 

-112.850878 

-112.795132 

-112.974426 

-112.979312 

-112.981302 

-112.7661185 

-112.946071 

-112.950276 

-112.952574 

2 

c 

gs 

.973 

.912 

.971 

.907 

.967 

.901 

R( 01 8-C4A) 
(bohrs ) 

4.6 

10.0 

SCF 

Cl 

EX 

DAV 

-112.703459 

-112.875073 

-112.878483 

-112.880629 

-112.65156 

-112.817392 

-112.819742 

-112.821163 

c" 

gs 

.966 

.902 

.973 

.905 

2'.  Oxetane  (OXET)  +  protonated  3-azidomethyl- 
3-methy loxetane  (AMM0H+) 

a’.  Results 

The  stabilization  point  R(01B-C4A)  * 
2.9  bohrs  and  19°  (fully  open)  with  the  stabilization  energy  -0.02386  a.u. 
-14.97  kcal/mole. 

The  potential  energy  surfaces  and  the  reaction  potential  map  are 
presented  in  Figures  11-27  to  11-29.  The  detailed  tables  of  results  follow 
in  Tables  11-34  to  11-38. 

Figure  11-27:  "MRO-CI  Extrapolated  Energy  for  OXET  Approaching  Protonated 
AMMO  For  Fixed  Angle  <5  and  Different  Intermodular  Distances 
R ( 01 B-C4 A ) " 

Figure  11-28:  "MRD-CI  Extrapolated  Energy  for  OXET-Protonated  AMMO  Complex 
For  Fixed  Intermol ecul ar  Distances  R(01B-C4A)  and  Different  6  Values" 
Figure  11-29:  "The  Potential  Energy  Surface  For  OXET  Approaching 
Protonated  AMMO" 

Table  11-34:  "OXET  +  AMM0H+  6=0°  (fully  closed)  Energies  (a.u.)  as  a 
function  of  R(01B-C4A)" 

Table  11-35:  "OXET  +  4MM0H+  <5-5°  Energies  (a.u.)  as  a  function  of 
R(018-C4A) " 

Table  11-36:  "OXET  +  AMM0H+  <5  -  10°  Energies  (a.u.)  as  a  function  of 
R ( 01 B— C4 A ) " 

Table  11-37:  "OXET  +  AMM0H+  5  -  15°  Energies  (a.u.)  as  a  function  of 
R ( 01 3-C4A) " 

Table  11-38:  "OXET  +  AMMOH'*’  <5  =  19°  (fully  open)  Energies  (a.u.)  as  a 
function  of  R(01B-C4A)" 


extrapolated  energy  for  oxet  approa 


R(0|B-  C4A) 


R ( 01 B-C4A) 
(bohrs ) 

2.3 

OXET  +  AMMOH 

6-0° 

(fully  closed) 

ENERGIES  (a.u.) 

2.9 

3.6 

SCF 

-112.164132 

-112.542090 

-112.706134 

Cl 

-L12. 358522 

-112.732924 

-112.890654 

EX 

-112.366976 

-112.739502 

-112.896460 

DAV 

-112.370250 

-112.742435 

-112.900469 

c2 

.964 

.9651 

.9569 

gs 

.905 

.906 

.906 

R(  01 B-C4A) 

(bo'nrs) 

4.6 

10.0 

SCF 

-112.757479 

-112.751545 

Cl 

-112.938395 

-112.933136 

EX 

-112.943841 

-112.936887 

OAV 

-112.945962 

-112.938947 

R ( 01B-C4A) 
(bohrs ) 

2.3 

OXET  +  AMMOH 

6  -  5° 

ENERGIES  (a.u.) 

2.9 

3.6 

SCF 

-112.380825 

—112. 654044 

-112.740730 

Cl 

-112.567865 

-112.840616 

-112.922895 

EX 

-112.574772 

-112.848736 

-112.928602 

DAV 

-112.577225 

-112.851497 

-112.932235 

c2 

.9687 

.965 

.958 

gs 

.9067 

.9057 

.9048 

R( 01 B-C4A) 

(bohrs ) 

4.6 

10.0 

SCF 

-112.755627 

-112.739371 

Cl 

-112.936087 

-112.921221 

EX 

-112.940998 

-112.924839 

DAV 

-112.943458 

-112.927221 

TABLE  11-36 


87 


% 

> 

V 


* 

1 


N 


c 

OXET  +  AMM0H+ 

6  -  10° 

ENERGIES  (a.u.) 


R ( 01 B— C4A) 
(bohrs) 

2.3 

2.9 

3.6 

SCF 

-112.517118 

-112.72324 

-112.75546 

Cl 

-112.697121 

-112.905683 

-112.927808 

EX 

-112.702417 

-112.911346 

-112.930962 

DAV 

-112.704430 

-112.913639 

-112.942933 

2 

c 

.909 

.9060 

.9124 

gs 

.971 

.9685 

.9226 

R( 01 B-C4A) 
(bohrs) 

4.6 

10.0 

SCF 

-112.73736 

-112.70814 

Cl 

-112.91068 

-112.874988 

EX 

-112.915738 

-112.876492 

DAV 

-112.91810 

-112.891173 

c2 

.904 

.909 

GS 

.965 

.9108 

% 

s 

■» 


■S 


E 


OXET  +  AMMOH 


ENERGIES  (a.u.) 


R ( 01 B-C4A) 
(bohrs ) 

2.3 

2.9 

3.6 

SCF 

-112.599563 

-112.763765 

-112.7622807 

Cl 

-112.774975 

-112.943217 

-112.940000 

EX 

-112.779670 

-112.948373 

-112.944830 

OAV 

C2 

-112.781491 

-112.950410 

-112.947124 

.973 

.970 

.9668 

GS 

.9118 

.9071 

.903 

R(0I B-C4A) 

(bohrs) 

4.6 

10.0 

SCF 

-112.7198416 

-112.677837 

Cl 

-112.889931 

-112.845810 

EX 

-112.894115 

-112.849130 

DAV 

-112.896185 

-112.850898 

*' 


OXET  +  AMMOH 

6  »  19° 
(fully  open) 

ENERGIES  (a. a.) 


R ( 01 B— C4 A) 
(bohrs ) 

2.3 

2.9 

3.6 

SCF 

-112.6397454 

-112.77667165 

-112.7547122 

Cl 

-112.811422 

-112.953525 

-112.932275 

EX 

-112.818284 

-112.960747 

-112.938784 

DAV 

-112.820049 

-112.962721 

-112.941037 

2 

r 

0.97329334 

0.97090597 

0.96709588 

v.. 

gs 

0.913082 

0.90737 

0.9018638 

Rf 01 B-C4A)  ...  . 

(bohrs) 

4.6 

10. U 

SCF 

Cl 

-112.866056 

-112.6445043 

-112.810255 

EX 

-112.871591 

-112.813953 

DAV 

-112.873703 

-112.815439 

c2 

0.96659015 

0.97269997 

as 

0.90336874 

0.90455537 

90 


3'.  3-Azidomethyl-3-methyloxetane  (AMMO)  + 
proton  a ted  3-azidomethy 1 -3-methy loxetane 

a'.  Results 


To  check  that  the  optimum  a  angle  for 

AMMO  +  AMM0H+  was  90°  we  ran  the  SCF  calculations  as  a  function  of  a. 

Table  11-39,  "AMMO  +  AMM0H+,  SCF  Energies,  a  Angle  Dependence"  The 
stabilization  point  R(C4A-01B)  -  2.9  bohrs  and  6  -  19°  (fully  open)  with 
stabilization  energy  =  -0.02408  a.u.  =>  -15.11  kcal/mol.  The  estimated 
activation  energy  is  6.27  kcal/mol. 

The  potential  energy  surfaces  and  the  reaction  potential  map  are 
presented  in  Figures  11-30  to  11-33.  The  detailed  tables  of  results  follow 
in  Tables  11-40  -  11-43. 


Figure  11-30:  "Extrapolated  Energy  For  AMMO  Approaching  Protonated  AMMO  For 
Fixed  Angle  6  And  Different  Intermolecular  Distances  R(01BC4A)" 

Figure  11-31:  "MR0-CI  Extrapolated  Energy  For  AMMO-AMMO  Protonated  Complex 
For  Fixed  Intermolecular  Distances  R(01B-C4A)  And  Different  6  Angle 
Values" 

Figure  11-32:  "AMM0-AMM0H*  Extrapolated  Cl  Energy  Along  the  Reaction  Coor¬ 
dinate  For  AMMO  Protonated  AMMO  Addition  Reaction" 


Figure  11-33:  "The  Potential  Energy  Surface  For  AMMO  Approaching  Protonated 
AMMO" 


Table  11-40:  "AMMO  +  AMM0H  +  ,  <5 
function  of  R(01B-C4A)" 

Table  11-41:  "AMMO  +  AMM0H  +  ,  <5 
R ( 01 B-C4A) " 

Tab1  e  11-42:  "AMMO  +  AMM0H  +  ,  <5 
3  ( 01 B — C  4  A ) " 

Table  1 1-43:  "AMMO  ♦  AMM0H  +  ,  <5 
R ( 01 B-C4A) " 

Table  11-44:  "AMMO  +  AMM0H+,  6 
function  of  R(018-C4A)" 


0°  (fully  closed)  Energies  (a.u.)  as  a 
5°  Energies  (a.u.)  as  a  function  of 
10°  Energies  (a.u.)  as  a  function  of 
15°  Energies  (a.u.)  as  a  function  of 
19°  (fully  open)  Energies  (a.u.)  as  a 


TRAPOLATED  ENERGY  FOR  AMMO  APPROACHING  PROTONATED  AMMO 
ED  ANGLE  ft  (A»0#,  B»6*,  C  =  10*,  D»15*f  E  *1 9° )  AND  DIFFERENT  INTERMOLECUL  AR 


R(0IB-C4A)(bohr) 


REACTION  COORDINATE 


R(OIB-C^A)  [bohr] 


TABLE  11-39 


AMMO  +  AMMOH+ 


a  ANGLE  DEPENDENCE 


a  (°) 


SCF 


0 

90 

130 

270 

305 

315 


-153.88041 

-153.97455 

-153.87978 

-153.97535 

-153.93036 

-153.93097 


TABLE  1 1 -a: 


AMMO  ♦  AMMOH 

6-0°  (FULLY  CLOSED) 
ENERGIES  (a.u.) 


R( 01 B-C4A) 
(bohrs ) 

2.3 

2.9 

3.6 

SCF 

-153.596296 

-153.975355 

-154.139594 

Cl 

-153.790525 

-154.165945 

-154.324062 

EX 

-153.798841 

-154.172725 

-154.330693 

DAV 

-153.802098 

-154.175653 

-154.333130 

c2 

0.964 

0.965 

0.967 

gs 

0.906 

0.906 

0.906 

R( 01 B-C4A) 
(bobrs) 

4.6 

10.0 

SCF 

-154.191083 

-154.186025 

Cl 

-154.371742 

-154.367499 

EX 

-154.377200 

-154.371299 

DAV 

-154.379303 

-154.373342 

c2 

0.970 

qs 

0.904 

TABLE  11-41 


AMMO  +  AMMOH* 

6  »  5° 

ENERGIES  (a.u.) 


R( 01 B-C4A) 
(bohrs ) 

2.3 

SCF 

-153.814107 

Cl 

-154.000988 

EX 

-154.007818 

OAV 

-154.010260 

c2 

0.968 

gs 

0.907 

R( 01 B-C4A) 

(bohrs ) 

4.6 

SCF 

-154.189225 

Cl 

-154.369528 

EX 

-154.374182 

OAV 

-154.376621 

2.9 

3.6 

-153.951247 

-154.174310 

-154.274199 

-154.356484 

-154.282183 

-154.362352 

-154.284928 

-154.365015 

0.966 

0.965 

0.906 

0.905 

10.0 

-154.173834 

-154.355578 

-154.359227 

-154.36160 

AMMO  +  AMMOH 


6  =  10u 

ENERGIES  (a.u.) 


R (01 3-C4A) 
(bohrs) 

2.3 

2.9 

3.6 

SCF 

-153.951248 

-154.157459 

-154.189200 

Cl 

-154. 130896 

-154.339691 

-154.367479 

EX 

-154.136065 

-154.345329 

-154.373005 

DAV 

-154.138029 

-154.347606 

-154.375491 

c2 

0.972 

0.969 

0.965 

gs 

0.809 

0.906 

0.904 

R( 01 B-C4A) 

(bohrs) 

4.6 

10.0 

SCF 

-154.170941 

-154.142564 

Cl 

-154.344119 

-154.316851 

EX 

-154.349152 

-154.320735 

DAV 

-154.351513 

-154.323112 

2 

c 

0.965 

0.963 

gs 

0.904 

0.900 

%.VA  JV-N  A  .%.v, 


S558® 


TABLE  11-43 


AMMO  +  AMMOH+ 

6  -  15° 

ENERGIES  (a.u.) 


R  ( 0 1 B— C4 A) 
(bohrs ) 

2.3 

SCF 

-154.034207 

Cl 

-154.209377 

EX 

-154.214183 

DAV 

-154.215995 

c2 

0.973 

gs 

0.912 

R( 01 B-C4A) 

(bohrs ) 

4.6 

SCF 

-154.153651 

Cl 

-154.323403 

EX 

-154.327627 

OAV 

-154.329684 

c2 

0.967 

gs 

0.905 

2.9 

3.6 

-154.198300 

-154.196284 

-154.376871 

-154.374171 

-154.381482 

-154.379496 

-154.384939 

-154.381804 

0.959 

0.967 

0.908 

0.903 

10.0 

-154.112237 

-154.280050 

-154.283300 

-154.285057 

0.968 

0.903 

TABLE  1 1 -4/1 


AMMO  ♦  AMMOH* 

5  »  19°( FULLY  OPEN) 
ENERGIES  (a.u.) 


R( 01 B-C4A) 
(bohrs ) 

2.3 

2.9 

3.6 

SCF 

-154.074638 

-154.211389 

-154.188879 

Cl 

-154.248331 

-154.390422 

-154.368442 

EX 

-154.253127 

-154.395379 

-154.373253 

DAV 

-154.254876 

-154.337350 

-154.375538 

c2 

0.973 

0.971 

0.967 

gs 

0.913 

0.907 

0.901 

R(013-C4A) 
(bohrs ) 

4.6 

10.0 

SCF 

-154.130128 

-154.078893 

Cl 

-154.301398 

-154.245371 

EX 

-154.305233 

-154.248153 

DAV 

-154.307349 

-154.249629 

c2 

0.966 

0.973 

gs 

0.903 

0.904 

(b).  Recap  of  Reaction  Energies  for  Cationic  Polymerization 
of  Energetic  Oxetanes  for  Initiation  and  Reaction 


Cationic  polymerization  has  two  major  steps: 
initiation  and  propagation.  Initiation  is  governed  by  the  propensity  for 
protonation  of  the  oxetane.  The  three  dimensional  electrostatic  molecular 
potential  contour  (EMPC)  maps  we  calculated  earlier  are  very  indicative  of 
the  propensity  of  the  energetic  substituted  oxetanes  to  initiate.  These 
EMPC  maps  are  also  indicative  of  the  propensity  of  the  energetic  substituted 
oxetanes  to  polymerize.  For  a  more  quantitative  comparison  of  propensity  to 
initiate  we  calculated  the  MRD-CI  energies  of  protonation  [AE(protonation) ] 
for  all  the  energetic  substituted  oxetanes  we  have  studied. 

The  next  step  in  cationic  polymerization  is  reaction  between  the 
oxetane  (or  substituted  oxetanes)  and  the  protonated  oxetane  (or  protonated 
substituted  oxetane).  We  have  calculated  the  MRD-CI  stabilization  energy 
[ AE ( add i tion) ]  for  several  series  of  reactants  as  a  function  of  the  angle 
(a)  between  the  rings,  the  interring  distance  R(01B-C4A)  and  the  angle  (6) 
of  opening  the  protonated  ring.  The  stabilization  point  for  all  of  the 
pairs  of  reactants  we  have  studied  to  date  is  R(01B-C4A)  -  2.9  bohrs  and  6  = 
19°.  (See  sketch  page  8  for  definition  of  6  angle) 

In  the  Table  (11-45)  are  tabulated  the  MRD-CI  values  for  AE 
(protonation),  AE  (addition)  and  AE  [AE  (protonation)  +  AE  (addition)]  at 

the  stabilization  point  for  OXET  (oxetane)  +  0XETH+  (oxetane  H+),  FNOX  + 

0XETH+,  OXET  +  FN0XH+,  FNOX  +  FN0XH+,  AMMO  +  0XETH+,  OXET  +  AMM0H+ ,  and  AMMO 

*  ammoh\ 

The  calculated  MRD-CI  values  for  AE  (protonation)  indicate  that  FNOX 
gains  less  energy  on  protonation  than  OXET  or  AMMO.  OXET  and  AMMO  have 
close  to  the  same  protonation  energy.  However,  our  earl ier  EMPC  maps  which 
indicated  that  OXET  had  a  larger  volume  within  its  -20  kcal  isopotential 
contour  would  imply  that  oxetane  would  still  have  a  somewhat  greater 
tendency  to  protonate  sooner  (initiate  sooner)  than  AMMO  even  though  they 
both  na/e  similar  AE  (protonation). 

The  calculated  AE  values  indicate  that  all  combinations  involving  FNOX 
or  f:N0XH+  are  less  favorable  than  those  involving  OXET  or  0XETH+  and  AMMO  or 
AMM0H+ . 


The  most  favorable  reactions  involve  0XETH+  reacting  with  OXET  or  AMMO 
rather  than  AMM0H+  reacting  with  OXET  or  AMMO. 

Our  initial  MRD-CI  result  on  protonation  of  BAMO  indicates  it  has  a 
lower  AE  (protonation)  than  that  of  AMMO.  Thus  in  a  BAM0-AMM0  mixture  AMMO 

will  tend  to  initiate  first  to  form  AMM0H+.  We  are  continuing  MRD-CI 


calculations  for  the  potential  surfaces  of  reactions  involving  BAMO  and 
3AM0H*  with  various  partners. 

The  results  of  calculations  such  as  these  will  enable  one  both  to 
understand  and  then  to  predict  copolymerization  preferences. 
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TABLE  11-45 

CATIONIC  POLYMERIZATION  INITIATION  AND  PROPAGATION 

PROTONATED  OXETANES  (OXETANEH*)  ♦  OXETANES 
AB-INITIO  MODPOT /VRDDO  MRD-CI 

ENERGIES  (a.u.) 


AE(protonation) 

AE(addi tion) 

AE 

OXETANE  + 

OXETANEH+ 

-0.31601 

-0.04378 

-0.35979 

FNOX  + 

OXETANEH* 

-0.31601 

-0.01113 

-0.32714 

OXETANE  + 

FNOXH+ 

-0.27068 

-0.06327 

-0.33394 

FNOX  + 

FNOXH+ 

-0.27068 

-0.03157 

-0.30225 

AMMO  + 

OXETANEH* 

-0.31601 

-0.04362 

-0.35963 

AMMO  + 

AMMOH+ 

-0.31548 

-0.02408 

-0.33956 

OXETANE  + 

AMMOH+ 

-0.31548 

-0.02386 

-0.33934 

OXETANE  * 

-0.30906 

8AM0H+ 


(c).  Population  Analyses 


Gerry  Manser  had  expressed  interest  in  the  population 
analyses  for  the  gross  atomic  populations  on  the  carbon  atom  being  attacked 
(C4A)  in  the  protonated  ring  and  the  total  overlap  populations  between  C4A- 
01 A  (in  the  protonated  ring)  and  01B-C4A  (the  bond  being  formed  between  the 
unprotonated  ring  and  the  protonated  ring). 

Following  are  detailed  tables  of  this  and  the  total  overlap  populations 

and  gross  atomic  populations  for  the  entire  AMMO  +  0XETH+,  OXET  +  AMM0H+, 

and  AMMO  +  AMM0H+  series  as  a  function  of  the  inter-ring  distance  R(01B-C4A) 
and  the  angle  of  opening  the  protonated  ring  (6).  Tables  11-46  -  Table  11-63 

The  general  conclusions  are: 

1'.  C2A  and  C4A  (the  a  carbons  in  the  original  protonated 
ring)  still  carry  about  the  same  excess  negative  charge 
(-0.2  e)  in  spite  of  the  fact  that  the  protonated 
species  carries  a  formal  positive  charge.  The  closeness 

of  the  charges  on  C2A  and  C4A  in  0XETH+  and  AMM0H+  is 
commensurate  with  their  similarities  in  pertinent  energy 
quantities.  As  the  protonated  ring  opens  and  the  01B  of 
the  unprotonated  ring  begins  to  form  a  bond  with  C4A, 

there  is  little  change  in  the  charge  on  C4A. 

2'.  01A  in  the  protonated  ring  (OXETH*  or  AMM0H+)  carries  an 

excess  negative  charge  of  -0.37.  01B  in  the 

unprotonated  ring  (OXET  or  AMMO)  carries  an  excess 
negative  charge  of  -0.35.  When  the  01B-C4A  bond  forms, 
018  still  carries  an  excess  negative  charge  of  -0.29  in 

AMMO  +  0XETH+  and  excess  negative  charge  of  0.32  in  AMMO 

+  AMM0H+.  When  the  C4A-01A  breaks,  the  01A  carries  a 

little  more  negative  charge  (0.44  in  AMMO  +  0XETH+,  0.45 

in  AMMO  +  AMM0H+)  in  both  protonated  partners. 

3'.  Total  overlap  populations  are  a  very  sensitive  criteria 
of  the  incipient  making  and  breaking  of  bonds.  The 
largest  TOP's  occur  when  the  energy  is  a  minimum. 

As  the  protonated  and  unprotonated  oxetane  (energetic 
substituted  oxetane) rings  approach 

a',  the  interring  TOP  (C4A-01A)  begins  to  get 

smaller  even  when  the  protonated  ring  is  still 
fully  closed.  This  indicates  that  the  C4A-01A 
bond  wants  to  lengthen. 


b‘.  The  TOP  (01B-C4A)  begins  to  be  noticeable  at 
4.6  bohrs  and  gets  larger  as  the  rings 
approach  closer  provided  that  the  protonated 
ring  is  open  by  at  least  6=5°.  The 
strongest  TOP ( 01 B-C4 A )  occurs  (as  anticipated) 
at  the  most  stable  point  energetically  R ( 0 1 3- 
C4A)  -  2.9  bohrs  and  6  -  19°. 
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8.  Ab-Initio  MRD-CI  Calculations  of  the  Protonation  of  Oxetane 

Protonation  of  oxetane  (or  substituted  oxetane)  is  the  initiation 
step  in  cationic  polymerization. 

We  had  long  since  pointed  out  the  initiation  protonation  step  in 
cationic  polymerization  is  an  ion-molecule  reaction 

A+  +  B  — >  AB+  (Type  II) 

where  the  ionization  potential  of  A  (the  hydrogen  atom  13.6  eV)  is  higher 
than  that  of  B  (oxetane  or  any  substituted  oxetane)  For  ion-molecule 
reactions  of  Type  I 

A+  ♦  B  — >  AB+  (Type  I) 

where  the  ionization  potential  of  A  is  lower  than  that  of  B  the  reaction  has 
the  possibility  to  proceed  along  a  single  lowest  potential  energy  surface 

which  separates  properly  theoretically  at  the  asymptotes  of  A+  and  B  (since 

both  A+  and  B  are  closd  shell  ground  state  systems).  However,  for  ion- 
molecule  reactions  where  the  ionization  potential  of  the  molecule  is  less 
than  13.6  eV  which  is  the  case  for  all  oxetanes  and  apparently  all  organic 
molecules)  there  has  to  be  at  least  one  potential  energy  surface  arising 

from  the  asymptotes  A  +  B+  which  is  lower  in  energy  at  the  asympotote  and 
along  at  least  part  of  the  interaction  potential  energy  surfaces. 

A  single  determinant  SCF  calculation  will  definitely  not  separate 
properly  at  the  asymptotes  and  is  not  sufficient  to  describe  a  protonation 
reaction  of  oxetanes  or  other  organic  molecules.  Ab-initio  MRD-CI 
calculations  for  the  entire  potential  energy  surface  for  ground  and  at  least 
several  excited  states  are  necessary  to  describe  protonation  reactions.  We 
had  long  stated  this  based  on  the  fundamental  physics  involved.  This  past 
year  we  verified  this  convincingly,  by  accurate  full  valence  electron  ab- 
initio  MRD-CI  calculations  for  the  ground  and  electronically  excitd  states 
for  the  proton  attack  on  oxetane.  We  carried  out  the  calculations  for  the 

linear  attack  of  the  proton  to  form  an  0-H+  bond  both  for  the  in-plane 
attack  and  the  out-of-plane  attack. 
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PROTONATION  OF  OXETANE 


Our  MRO-CI  results  confirmed  our  earlier  geometry  optimization  results 

that  the  0-H+  bond  is  in  the  plane  of  the  oxetane  ring  even  though  the  lone 
pairs  on  the  original  oxetane  are  out  of  plane.  The  results  indicated  an 
even  greater  complexity  for  the  potential  energy  surfaces  then  we  had 
originally  hypothesized.  The  lowest  energy  root  from  1.5  bohrs  to  4.5  bohrs 

(including  the  equilibrium  internuclear  0-H+  distance  of  protonated  oxetane 

at  2.0  bohrs)  is  *A,.  At  4.5  bohrs  the  lowest  root  becomes  to  an 

+  2  1 
asymptote  of  oxetane  (ground  state,  B^)  +  H.  At  4.5  bohrs  the  A^  curve 

continues  as  the  second  root  of  the  Cl  matrix  smoothly  to  oxetane*(first 

2  1 
excited  singlet  state  A^)  ♦  H.  The  third  root  of  the  Cl  matrix,  also  a  A1 

state,  has  a  minimum  at  2.0  bohrs,  a  hump  at  2.75  bohrs  and  then  continues 

down  to  the  separated  products  oxetane(ground  state)  +  H+.  The  behavior  of 
the  next  two  higher  roots  of  the  Cl  matrix  is  also  complicated.  This 

implies  that  while  the  ^  potential  energy  surface  arising  from  or 

dissociating  to  oxetane  (or  substituted  oxetane)  +  H+  will  always  be  one  of 
two  higher  lying  potential  energy  surfaces  involving  protonated  oxetane  (or 
protonated  substituted  oxetane),  the  exact  relative  position  of  that 
potential  energy  surface  will  be  critically  dependent  on  the  explicit  values 


of  the  ionization  potentials  for  the  substituted  oxetane  (ground  or  excited 
states)  and  will  also  depend  on  the  type  of  molecular  orbital  from  which 
ionization  had  occurred  (i.e.  lone  pair,  bonding  or  non-bonding  orbital). 
Thus  the  behavior  upon  protonation  can  be  somewhat  different  for  each 
individual  substituted  oxetane. 

These  findings  have  a  profound  significance.  Thus  MRD-CI  calculations 
for  the  ground  and  a  number  of  excited  states  will  be  necessary  to  describe 
protonation  or  deprotonation  processes. 

This  behavior  will  hold  true  for  all  protonation  reaction  of  organic 
molecules  since  their  ionization  potentials  are  less  than  13.6  eV.)  as  well 
as  protonation  reactions  involving  molecules  with  ionization  potentials 
lower  than  the  ionization  potential  of  the  hydrogen  atom,  13.6  eV.  This 
multipotential  surface  behavior  also  holds  true  for  deprotonation  of  these 
protonated  species. 

Thus  the  potential  energy  surfaces  for  protonation  reactions  are  much 
more  complicated  than  customarily  assumed  and  moreover  are  crucially 
dependent  not  only  on  the  lowest  ionzation  potentials  of  the  molecules 
involved  but  also  on  the  energies  of  their  higher  ionization  potentials 
relative  to  the  ionization  potential  of  the  hydrogen  atom. 


III.  Ab-Initio  MRD-CI  Calculations  for  Breaking  a  Chemical  Bond  in  a 

Molecule  in  a  Crystal  or  Other  Solid  Environment 

Breaking  a  > C-NO2  or  >N-N02  bond  is  the  initial  step  in  fractoemi ssion 

of  explosives  and  also  the  initial  step  leading  to  detonation  of  explosives. 
To  describe  properly  breaking  of  a  chemical  bond  in  a  molecule  it  is 
necessary  to  to  carry  out  ab-initio  MRD-CI  (mul tireference  double  excitation 
configuration  interaction)  calculations  of  the  isolated  molecule.  To 
describe  properly  breaking  of  a  chemical  bond  in  a  molecule  in  a  crystal  or 
other  solid  environment  it  is  necessary  to  carry  out  ab-initio  MRD-CI 
calculations  on  the  molecule  surrounded  by  other  molecules  as  in  the  crystal 
or  solid  arrangement.  Even  this  generation  of  supercomputers  still  does  not 
have  the  space  to  carry  out  such  calculations  on  large  nitroexplosive 
molecules  especially  since  many  of  them  (such  as  RDX  and  HMX)  have  a  large 
number  of  molecules  in  the  unit  cell. 

We  had  earlier  derived  and  implemented  and  used  successfully  a  method 
for  dissociation  of  large  molecules  based  on  localized/local  orbitals.  The 
localized  molecular  orbitals.  The  localized  molecular  orbitals  in  the 
region  of  the  bond  breaking  are  included  explicitly  in  the  MRD-CI.  The 
remainder  of  the  occupied  and  virtual  orbitals  are  folded  into  an 
"effective"  Cl  Hamiltonian. 

A.  Methodology 

MRO-CI  calculations  are  absolutely  necessary  to  describe 
bondbreaking  processes  correctly  in  the  ground  state  and  especially  in  the 
excited  states. 

The  electronic  challenge  arose  to  extend  our  MRD-CI  (multireference 
double  excitation  -  configuration  interaction)  technique  based  on 
localized/local  orbitals  to  the  breaking  of  a  chemical  bond  in  a  molecule  in 
crystal  (or  other  solid  environment).  This  past  year  we  have  derived, 
implemented,  and  used  successfully  a  procedure  for  doing  this.  We  made  the 
first  presentation  of  results  using  this  method  spring  1987  at  the  ONR 
Workshop  on  Dynamic  Deformation,  Fracture  and  Transient  Combustion  of 
Energetic  Compounds. 

Our  technique  involves  solving  a  quantum  chemical  ab-initio  SCF 
explicitly  for  a  system  of  a  molecule  surrounded  by  a  number  of  other 
molecules  (the  unit  reference  cell  or  larger  assemblage)  in  the  multipole 
environment  of  yet  more  further  out  surrounding  molecules.  Multipoles  in 
the  environmental  region  affect  the  one-electron  term  in  the  Hamiltonian. 
This  Hamiltonian  is  solved  for  the  SCF  for  all  the  molecules  in  the  space 
treated  explicitly  quantum  chemically.  The  resulting  canonical  molecular 
orbitals  are  localized.  All  of  the  occupied  and  virtual  localized  orbitals 
in  the  region  of  interest  are  included  explicitly  in  the  MRD-CI  and  the 
remaining  occupied  localized  orbitals  are  folded  into  an  "effective"  Cl 
Hamiltonian.  The  advantage  is  that  the  transformations  from  integrals  over 
atomic  orbitals  to  integrals  over  molecular  orbitals  (the  computer  time-, 
computer  core-  and  external  storage  -  consuming  part  of  the  Cl  calculations) 
only  have  to  be  carried  out  for  the  localized  molecular  orbitals  included 
explicitly  in  the  MRD-CI  calculations. 
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Space  is  broken  up  into  three  regions: 

(C[B[A]B']C) 

A  Localized  space  treated  explicitly  in  ab-initio  MRD-CI 

calculations.  (This  can  be  an  entire  molecule  or  the  localized 
dissociation  region  of  a  large  molecule.) 

B  +  A  +  8'  Space  treated  explicitly  quantum  chemically  (ab-initio  SCF)  for 
supermolecule  BAB' 

C  +  C‘  Space  represented  by  multipoles  of  additional  molecules  taken 
into  account  by  inclusion  of  multipole  interactions  (up  through 
quadruooles)  into  one-electron  part  of  SCF  Hamiltonian. 

This  method  is  completely  general.  The  space  treated  explicitly 
quantum  chemically  and  the  surrounding  space  can  have  defects,  deformations, 
dislocations,  impurities,  dopants,  edges  and  surfaces,  boundaries,  etc. 

To  be  able  to  carry  out  such  MRD-CI  calculations  for  breaking  a 
chemical  bond  in  a  molecule  or  a  crystal  (or  other  solid  environment) 
represents  a  significant  breakthrough. 

The  desirable  optimal  computational  strategies  we  have  developed  over 
the  years  for  ab-initio  calculations  on  large  molecules  and  molecular 
systems  is  what  makes  these  computations  tractable  for  us:  MODPOT  -  ab- 
initio  effective  core  model  potentials  which  enable  calculation  of  the 
valence  electrons  only  explicitly,  yet  accurately;  VRDDO  -  a  charge- 
conserving  integral  prescreening  evaluation,  especially  effective  for 
spatially  extended  systems;  MERGE  -  to  save  and  reuse  integrals  from  common 
fragments  or  molecules,  which  enables  us  to  build  up  larger  and  larger 
clusters  very  efficiently  and  a  special  SCF  technique  which  allows  us  to  use 
the  SCF  wave  functions  from  smaller  clusters  as  a  start  to  get  rapid 
convergence  for  larger  clusters. 

Test  results  on  5  unit  cells  of  H£  (for  dissociating  the  molecule  in 

the  center  cell  treated  by  this  method)  verified  the  validity  of  this  new 
technique. 

B.  Calculations  Carried  Out  for  Nitromethanes : 


We  have  been  carrying  out  extensive  test  calculations  by  this  new 
technique  for  the  dissociation  of  the  H-jC  -  NO2  bond  in  nitromethane  for 

various  numbers  of  molecules  treated  explicitly  on  the  SCF  in  the  multipole 
field  of  varying  numbers  of  additional  CH^NO^  molecules  as  in  the  crystal 

arrangement  followed  by  localization  and  ab-initio  MRD-CI  calculations  on 
breaking  the  CH^  -  NO2  bond  in  a  specific  nitromethane  molecule.  Since  this 


m 


technique  is  new  we  are  carrying  out  extensive  testing  to  ascertain  how  many 
molecules  must  be  treated  in  each  region  for  reliable  results. 

One  of  the  pertinent  questions  we  posed  initially  for  decomposition  of 
molecules  in  crystals  was  did  it  take  more  or  less  energy  to  break  the  bond 
when  the  molecule  was  in  a  crystal  compared  to  breaking  the  bond  of  an 
isolated  molecule.  The  MRD-CI  results  for  breaking  the  H^C  -  NO^  bond  of 

nitromethane  in  the  presence  of  multipoles  of  other  nitromethane  molecules 
compared  to  the  MRD-CI  results  for  breaking  the  H^C  -  NO,,  bond  in  an 

isolated  nitromethane  molecule  indicate  that  it  takes  more  energy  to  break 
the  H-jC  -  N0?  bond  when  nitromethane  is  in  the  field  of  the  additional 

nitromethane  molecules. 

We  have  also  investigated  such  MRD-CI  calculations  treating  different 
numbers  of  nitromethane  molecules  explicitly  in  the  SCF  calculation  and 
varying  the  numbers  of  the  external  nitromethane  molecules  contributing  to 
the  multipole  field. 

Several  different  choices  of  the  arrangement  of  nitromethane  molecules 
was  also  considered  picking  different  pieces  from  the  experimental  crystal 
structure  of  ni tromethane. 

The  ab-initio  MRD-CI  calculations  have  been  carried  out  for  Rr  .. 

C-N 

distances  for  the  H^C-NC^  molecule  being  dissociated  of  2.4,  2.8,  3.0,  3.2, 
3.6,  4.4  and  5.6  bohrs  for  each  choice  of  nitromethane  system  description. 


C.  Detailed  Results  of  Calculations  Carried  Out  for  Nitromethane: 

Various  Choices  of  Size  and  Description  of  System. 

All  Figures  and  Tables  appear  at  the  end  of  this  section  starting  on 
page  136 

1.  The  Model  Case 

To  test  the  proposed  approach  the  model  cluster  has  been 
selected  from  real  nitromethane  crystal.  The  model  cluster  includes: 

a.  reference  molecule  (A) 

b.  two  closest  neighbors  (B)  to  reference  molecule  (these  two 
molecules  are  the  most  important  for  decomposition  of  C  -  N  bond  in  the 
reference  molecule) 

c.  two  far-lying  molecules  (C)  having  big  influence  on  the  B 
region  (Fig  III-l,  "The  Cluster  of  5  Nitromethanes  Chosen  For  Model 
Calculations) " . 

To  test  how  good  is  the  representation  of  a  molecule  by  multipoles, 
calculations  for  two  different  geometries  of  the  cluster  have  been 
performed.  (See  Figure  1 1 1 - 1 )  First  full  SCF  calculation  for  three 
nitromethane  molecules  with  the  reference  nitromethane  molecule  at  position 
A  and  two  more  nitromethane  molecules  at  positions  B  compared  to  SCF 
calculation  for  the  reference  nitromethane  molecule  at  position  A  and  two 
nitromethane  molecules  represented  by  multipoles  at  position  B.  Second, 
full  SCF  calculation  for  three  nitromethane  molecules  with  the  reference 
molecule  at  position  A  and  two  nitromethane  molecules  at  position  C  compared 
to  SCF  calculation  for  the  reference  nitromethane  molecule  at  position  A  and 
two  nitromethane  molecules  represented  by  multipoles  at  position  C.  These 
results  are  presented  in  Table  III-l,  "Energies  For  Equilibrium  Geometry 
(Rcn  -  3.0  bohrs),  Bond  Dissociation  Energy  and  Relative  Error  For  Two 

Geometries  For  the  Cluster  of  3  Nitromethane  Molecules".  The  results  show 
that  the  multipole  representation  is  extremely  reliable  compared  to  full  ab- 
initio  SCF  results  for  the  three  nitromethane  molecules.  The  relative  error 
on  the  MRD-CI  calculated  C-N  dissociation  energy  when  the  two  nitromethane 
molecules  in  the  C  position  are  represented  by  their  multipole  effects  is 
only  0.027%.  The  relative  error  when  the  two  closer  nitromethane  molecules 
in  the  B  positions  are  represented  by  their  multipole  effects  is  still  only 
1.14*. 

Table  1 1 1 -2 ,  "Total  Energies,  Reduced  Energies,  and  Bond  Dissociation 
Energies  For  Model  Cluster  of  Nitromethane  Molecules  Within  Different 


Approaches",  has  more  extensive  results  for  the  E^p, 
dissociation  energy  for: 


-Cl, EX 


-N  bond 


5  n i tromethanes  fun  SCF 

3  n i tromethanes  full  SCF  in  field  of  2  multipoles 

3  ni tromethanes  full  SCF 

1  nitromethane  full  SCF  in  field  of  4  multipoles 

1  nitromethane  full  SCF  in  field  of  2  multipoles  in  B  positions 

1  nitromethane  single  free  molecule  full  SCF 


Table  1 1 1 -2a ,  "Total  Energies  and  Bond  Dissociation  Energies 


(Corresponding  to  R 


CN 


5.6  bohrs  and  R 


CN 


®)  For  Model  Cluster  of 


Nitromethane  Molecules  Within  Different  Approaches",  presents  results  for 

CN 

cluster  as  in  Table  1 1 1 -2 . 


R^  =*  5.6  bohrs  and  Rrw  =  ®  for  the  same  descriptions  of  the  nitromethane 


The  bond  dissociation  energy  of  a  molecule  in  a  crystal,  (or  other 
solid)  in  contrast  to  the  bond  dissociation  energy  in  a  free  molecule, 
includes  interactions  with  other  molecules  in  the  crystal  or  other  solid. 
The  energy  of  the  cluster  before  the  decomposition  can  be  written 


'1 


ea  +  eb+  eab+  eac+  ecb+  EABC’ 


( A, B,C  correspond  to  spaces  of  cluster,  Fig  1 1 1 - 1 .  where  here  B  space  =  B 
B* ; 

C  »  C  +  C‘)  and  after  decomposition 


E«  +  E  +  Er»A  + 
A  B  BC 


where  E^,  E^  energy  of  molecule  A  before  and  after  decomposition 


E.  . 
TJ 


Eijk 


energy  of  .molecules  B 

(Eg  -  Eg  +  Eg,  +  Egg,  in  our  case) 

two  body  interactions 

th  ee  body  interactions 


The  bond  dissociation  energy 


AE  ’  el  “  e2  ■  el  “  EA  •  EB 


Assuming  the  decomposition  of  the  bond  to  infinity,  E^  is  a  energy  of 
completely  decomposed  free  molecule  and  ER  is  an  energy  of  the  cluster 


without  A  molecule. 


The  value 

E^( r )  =  f )  ~  Eg, 

is  called  the  reduced  energy  and  represents  the  energy  of  reference  molecule 
in  the  field  of  other  spaces  of  the  crystal.  The  ER(r)  can  be  used  to 

compare  energy  surfaces  (reaction  surfaces)  in  different  assumed  models. 

These  calculations  have  also  been  carried  out  as  a  function  of  R~  ...  Table 

L-N 

III-3,  "The  Decomposition  Pathway  for  C-N  Bond  in  Nitromethane  Crystal 
(Model  Cluster)  For  Different  Approaches" 

2.  Real  Cluster  -  Extended  Cluster  and  Alternative  Extended  Cluster 

The  cluster  of  13  molecules  has  been  chosen  to  represent  the 
real  crystal.  Within  this  cluster  the  distance  between  atoms  C^  and  of 

reference  molecule  and  C  or  N  atom  in  neighbor  molecules  is  smaller  than  10 
bohrs.  Since  only  5  ni tromethanes  can  be  treated  quantum  mechanically  at 
present  due  to  computer  disc  limitations  (which  limits  the  size  of  the 
transformation  from  atomic  orbitals  to  molecular  orbitals  for  the  MRD-CI 
calculation) ,  the  selection  of  4  molecules  around  the  reference  molecule  was 
necessary.  Two  choices  have  been  tried.  First:  four  molecules  closest  to 
reference  molecule  (Fig.  III-2,  "The  Cluster  of  5  Nitromethanes  From 
Extended  Cluster")  and  second  four  molecules  having  atoms  which  are  the 
closest  to  C-N  bond  (Fig.  III-3,  "The  Cluster  of  5  Nitromethanes  From 
Extended  Cluster  (Alternative  Choice").  Eight  remaining  molecules  were 
represented  by  multipoles.  Multipoles  were  generated  in  crystal  orbital 
calculation),  by  the  crystal  orbital  part  (XTLORB)  of  our  P01Y-CRYST 
program  (for  ab-initio  calculations  on  polymers  and  crystals)  taking 
advantage  of  repeated  symmetry  units  for  one  unit  cell  and  then  decomposed 
into  point  charges. 


a.  Extended  Cluster 
Figure  1 1 1 - 2 

Table  1 1 1 -4 :  "Total  Energies,  Reduced  Energies,  and  Bond 
Dissociation  Energies  for  Extended  Cluster  of 
Nitromethane  Molecules  Within  Different  Approaches 
(energy  in  a.u.,  equilibrium  distance  R  -  3.0  bohrs)" 

Table  III -4a :  "Total  Energies,  Reduced  Energies,  and  Bond 
Dissociation  Energies  (corresponding  to  -  5.6 

bohrs  and  R^  -  ®)  For  Extended  Cluster  of 

Nitromethane  Molecules  Within  Different  Approaches" 


b.  Extended  Cluster  (Alternative  Choice) 

Figure  III -3 

Table  1 1 1 -5  "Total  Energies,  Reduced  Energies,  and  Bond 
Dissociation  Energies  for  Extended  Cluster  of 
Nitromethane  Molecules  (Alternative  Choice)  Within 
Different  Approaches  (energy  in  a.u.,  equilibrium 
distance  R  -  3.0  bohrs) 

Table  III -5a :  "Total  Energies,  Reduced  Energies,  and  Bond 
Dissociation  Energies  (corresponding  to  RCN  -  5.6 

bohrs  and  R^N  -  «)  For  Extended  Cluster  of 

Nitromethane  Molecules  (Alternative  Choice)  Within 
Different  Approaches" 

For  both  the  extended  cluster  and  the  extended  cluster  (alternative 
choice)  these  calculations  have  also  been  carried  out  as  a  function  of  RC_N. 

Table  1 1 1 -6 :  "The  Decomposition  Pathway  for  C-N  bond  in 

Nitromethane  Crystal  (Extended  Cluster)  for  Different 
Approaches." 

Table  1 1 1—7 :  "The  Decomposition  Pathway  for  C-N  bond  in 
Nitromethane  Crystal  (Extended  Cluster,  Alternative 
Choice)  For  Different  Approaches" 
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3.  Crystal  Orbital  Calculations  For  Model  Nitromethane  Cluster 

To  check  the  goodness  of  the  model  cluster  approximation  for 
crystalline  nitromethane  we  carried  out  ab-initio  crystal  orbital 
calculations  (XTLQRB)  using  our  POLY-CRYST  program.  The  POLY-CRYST 
program  calculates  ab-initio  crystal  orbitals  for  the  unit  reference  cell 
(or  unit  reference  polymer  segment)  taking  advantage  of  the  translational 
symmetry  in  a  crystal  and  the  t rans 1  at ional /rotational  symmetry  in  a 
polymer.  This  POLY-CRYST  program  incorporates  as  options  all  the  desirable 
computational  strategies  we  had  derived  over  the  years  for  calculations  on 
large  molecules:  Ab-initio  M00P0T,  VRDDO,  and  MERGE.  (See  description  page  128 
III.  Section  A.  Methodology  Ab-Initio  MODPOT/VRDDO/MERGE. )  These  XTLORB 
calculations  were  carried  out  for  three  cells  of  4  nitromethane 
molecules/unit  cell  (from  the  model  cluster).  These  XTLORB  calculations 

were  carried  out  integrating  over  3  £  points  and  over  5  £  points.  Results 
(Table  III -8 ,  "Nitromethane  -  Crystal  Orbital  Calculations  for  Model 
Cluster")  showed  the  XTLORB  calculations  for  nitromethane  integrating  over 

only  3  i<  points  were  essentially  already  converged.  This  demonstrates  that 

for  molecular  crystals  integration  over  far  fewer  £  points  is  necessary  for 
convergence  than  for  atomic,  interatomic  or  ionic  crystals.  These  XTLORB 
calculations  were  then  carried  out  for  3  cells  of  three  nitromethane 
molecules/unit  cell  (taking  one  central  molecule  away  from  the  cluster). 

The  difference  in  the  XTLORB  total  energies  between  the  4  nitromethane 
molecules/unit  cell  and  the  3  nitromethane  molecules/unit  cell 
(Table  III -8 . ) 

Er  -  E4  -  E3  -  -48.06091  a.u. 

corresponds  very  closely  to  the  reduced  energy  per  nitromethane  molecule, 

ER  »  -48.06057  a.u. 

calculated  from  explicit  SCF  calculations  on  the  model  nitromethane  cluster 
in  the  multipole  field  of  farther  out  nitromethane  molecules  for  the  model 
cluster. 

Thus,  the  multipole  approximation  for  describing  the  effect  of  further 
out  molecules  on  the  SCF  cluster  energies  is  quite  good. 


D.  Conclusions 

To  be  able  to  carry  out  such  MRD-CI  calculations  for  breaking  a 
chemical  bond  in  a  molecule  or  a  crystal  (or  other  solid  environment) 
represents  a  significant  breakthrough . 

Describing  the  processes  and  mechanisms  of  breaking  chemical  bonds  in 
an  energetic  molecule  when  the  molecule  is  in  a  crystal  or  other  solid 
environment  leads  to  an  understanding  of  the  initiation  of  energetic 
reactions  and  subsequent  processes  leading  to  detonation  and  also  leads  to 
an  understanding  of  fractoemi ssion  processes. 

The  results  confirm  the  preliminary  results  of  our  original  MRD-CI  test 
calculations  for  breaking  the  H^C  -  NO,,  bond  in  a  single  nitromethane  in  the 

multipole  field  of  6  other  nitromethane  molecules,  namely  that  it  takes  more 
energy  to  break  the  H^C  -  NO^  bond  when  the  nitromethane  molecule  is  in  a 

crystal  or  solid  environment  in  the  presence  of  other  nitromethane 
molecules.  However,  to  reach  the  convergent  value  it  is  necessary  to  take 
into  account  a  number  of  nitromethane  molecules  explicitly  in  the  SCF  and  a 
number  more  nitromethane  molecules  in  the  surrounding  multipole  space.  We 
are  continuing  our  investigations  in  this  area. 

Our  development  has  led  to  crucial  understanding  of  the  initial  steps 
leading  to  detonation.  We  plan  to  explore  the  application  of  our  techniques 
to  the  subsequent  events  taking  place  after  the  initial  bond  breaking. 

We  plan  to  carry  out  such  calculations  for  breaking  the  >C-N02  and  >N- 

NO2  bond  breaking  in  larger  energetic  molecules  such  as  RDX  and  other 

energetic  compounds.  This  method  is  very  applicable  for  such  calculations. 
Carrying  out  such  calculations  on  larger  molecules  will  be  dependent  on  the 
computer  core,  disc  storage  and  time  available.  However,  computers  are 
getting  faster  and  their  memories  and  peripheral  storage  are  getting  larger. 
The  next  generation  of  computers  which  is  already  past  the  drawing  board 
stage  are  ideally  suited  for  such  calculations. 
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the  cluster  OF  5  nitromethanes  from  extended 

CLUSTER  (ALTERNATIVE  CHOICE) 
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TABLE  III-l 


Energies  For  Equilibrium  Geometry  (R^  =•  3.0  bohrs).  Bond 

Dissociation  Energy  and  Relative  Error  For  Two  Geometries  For  The 
Cluster  of  3  Nitromethane  Molecules  (see  Fig  1 1 1  —  1 )  (energy  in 
[a.u.]) 

CLUSTER  Ej-£p  Eqj  C-N  bond  diss-  relative 

dissociation  energy,  error  [%] 
AE 


3  NITROMETHANES 

C  -  A  -  C' 

1  NITROMETHANE 
+  2  MULTIP 
in  position  C 


el  -144.136254 
1) 

Er  -48.051505 
-48.052639 


-144.388879 

-48.304130 

-48.304109 


.0077714 

.077683 


0.027% 


3  NITROMETHANES 

8  -  A  -  B* 

1  NITROMETHANE 
+  2  MULTIP 
in  position  B 


ej  -144.146063 
2) 

Er  -48.061789 
-48.062455 


-144.396239 

-48.311965 

-48.312941 


.085549 

.086525 


1)  Er  -  e,  -  Escp  (C  -  C') 
2>  er-£i  -escf  (b-b<) 


-1.14% 
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TABLE  1 1 1-2 


Total  Energies,  Reduced  Energies,  and  Bond  Dissociation  Energies 
for  Model  Cluster  of  Nitromethane  Molecules  Within  Different 
Approaches  (energy  in  [a.u.]) 


Approach 

escf 

ECI , EX  C_N  bond 

dissociation 
AE 

5  NITROMETHANES 

e1  -240.238726 

-240.489022 

1) 

ER  -48.060579 

-48.310875 

-0.084458 

3  NITROMETHANES 

e.  -144.154520 

-144.404582 

+  2  MULTIP 

2) 

Er  -48.060636 

-48.310698 

-0.084282 

3  NITROMETHANES 

e,  -144.146063 

-144.396239 

3) 

Er  -48.0617891 

-48.311965 

-0.085549 

NITROMETHANE 
+  4  MULTIP 

-48.061205 

-48.311927 

-0.085511 

NITROMETHANE 
+  2  MULTIP 
in  position  B 

-48.062455 

-48.312941 

-0.086525 

NITROMETHANE  -48.052639 

SINGLE  FREE  MOLECULE 

-48.305044 

-0.078628 

ER  "  el  escf 

(4  NITROMETHANES, 

CBB'C' ) 

2)  ER  “  el  escf 

(2  NITROMETHANES, 

BB ' 

♦  2  MULTIP  C,C' 

’) 

3)  ER  "  el  '  escf 

(2  NITROMETHANES, 

BB') 

W  l»  1"  w l’  u u  IfW  ».■  u 


1  *."  I."  r  V>.  V  V  nr  »rv  -v  v.  V"v  \  -v 


.*v.  *v  ■  m  j-w  --rw  v 
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TABLE  III -8 


Ni tronetnane  -  Crystal  Orbital  Calculations  for  Model  Cluster 

Crystal  Orbital  Calculations  for  4  Molecules  3  Cell  Calculations  in  i  Ax 
Oirection  (Fig.  1 1 1 - 1  ) 

Total  Energy  for  4  molecules 

3  k  points  E4  -  -192.195942921  a.u. 

5  ic  points  E4  -  -192.195942922  a.u. 

Total  Energy  for  3  molecules 

3  k  points  E^  «  -144.13503504  a.u. 

5  k  points  E4  -  -144.13503504  a.u. 

ER  -  E4  -  E3  -  -48.0609079  a.u. 

P  energy  from  cluster  calculations 
Er  -  -48.060579 


IV.  POL Y-CRYST 


POLY-CRYST  is  the  program  we  previously  derived  and  wrote  for  ab-initio 
calculations  on  crystals  and  polymers  using  the  translational  symmetry  in  a 
crystal  and  the  translational-rotational  symmetry  in  a  polymer.  Commensurate 
with  the  ONR  priorities  expressed  to  us  by  our  ONR  Contract  Monitor,  we 
devoted  only  minimal  but  still  scientifically  significant  effort  to  further 
development  and  testing  of  the  POLY-CRYST  program.  As  options  we  had 
already  included  POLY-CRYST  our  own  ab-initio  MODPOT  (ab-initio  effective 
core  model  potentials)  and  VRODO  (a  charge  conserving  integral  prescreening 
procedure)  options.  It  is  these  features  particularly  VRDDO  which  enables 
POLY-CRYST  to  handle  molecular  crystals  of  large  molecules  and  with  large 
numbers  of  large  molecules  per  unit  cell.  This  year  we  derived  and 
incorporated  into  POLY-CRYST  including  the  multipole  effects  of  farther  out 
molecules  to  take  into  account  long  range  effects  also.  We  then  meshed  this 
multipole  procedure  back  into  the  MRD-CI  programs  to  enable  us  to  include 
multipole  effects  when  breaking  a  chemical  bond  in  a  crystal.  We  also  ran 
some  tests  on  POLY-CRYST  on  integral  thresholds  and  numbers  of  unit  cells 
necessary  for  convergence.  These  preliminary  tests  identified  necessary 
criteria. 

Toward  this  convergence  criteria  goal,  we  also  derived  and  implemented 
a  procedure  for  calculating  the  charge  imbalance  caused  by  various  integral 
thresholds  to  give  a  precise  measure  of  the  effect  on  the  crystal  orbital 
calculation  of  dropping  integrals  of  various  sizes.  The  POLY-CRYST  program 
has  promise  for  yielding  important  fundamental  results  on  crystalline 
energetic  materials. 

As  mentioned  in  Section  III  we  also  used  this  XTLORB  program  to 
calculate  a  3  cell  case  of  nitromethane  to  verify  the  validity  and  accuracy 
of  the  procedure  we  use  for  the  effect  of  multipoles  from  farther  out 
molecules  on  SCF  calculations  for  molecular  clusters  from  crystals  the  SCF 
calculation  is  carried  out  explicitly  for  all  the  molecules  in  a  unit 
reference  cell  (or  larger  piece  of  the  crystal).  Our  SCF  results  for  the 
nitromethane  molecular  cluster  from  the  nitromethane  crystal  field  of  yet 
further  out  molecules  were  very  close  to  those  full  XTLORB  results  for 
nitromethane.  The  SCF  orbitals  are  then  localized  and  the  MRD-CI 
calculations  are  carried  out  for  breaking  a  chemical  bond  in  a  molecule  in  a 
crystal  including  explicitly  in  the  MRD-CI  calculation  the  localized 
orbitals  in  the  region  of  interest. 
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V.  Lectures  Presented  and  Publications  on  This  ONR  Research 


Presentations  given  and/or  scheduled  and  papers  published  and/or 
submitted  during  the  fiscal  year. 

A.  Presentations  Given  Or.  Joyce  J.  Kaufman 


Already  Presented 


(*  denotes  invited  lecture) 


a.  At  National  and  International  Meetings 

*  (1).  "Ab-Initio  MRD-CI  Calculations  for  the 
Propagation  Step  in  the  Cationic  Polymerization  of  Oxetanes  Based  on 
Localized  Orbitals,"  an  invited  paper  presented  at  the  International 
Symposium  on  Atomic,  Molecular  and  Solid  State  Theory,  Marineland,  Florida, 
March  1987. 

*  (2).  "Comparison  of  Ab-Initio  M00P0T  and  Ab-Initio 
Energy  Partitioned  Potential  Functions  for  Nitromethane  Dimer  Against  Large 
Basis  Set  Calculation,"  an  invited  paper  presented  at  the  Sanibel 
International  Symposium  on  Atomic  Molecular  and  Solid  State  Theory, 
Marineland,  Florida,  March  1987. 

*  (3).  "Ab-Initio  MRD-CI  Calculations  on  Protonated 
Cyclic  Ethers.  I.  Protonation  Pathways  Involve  Multipotential  Surfaces 
(Protonation  of  Oxetane)  II.  Differences  from  SCF  in  Dominant  Configurations 
Upon  Opening  Non-Protonated  Oxirane  Rings  (Epoxides),"  an  invited  paper 
presented  at  the  Sanibel  International  Symposium  on  Atomic  Molecular  and 
Solid  State  Theory,  Marineland,  Florida,  March  1987. 

*  (4).  "Ab-Initio  Calculations  on  Large  Molecules  and 
Solids  Using  Desirable  Computational  Strategies,"  Vlllth  International 
Conference  on  Computers  in  Chemical  Research  and  Education,  Beijing,  China, 
June  1987. 

(5).  "Ab-Initio  Quantum  Chemical  Calculations  on 
Large  Molecular  Systems  and  Crystals,"  American  Conference  on  Theoretical 
Chemistry,  Gull  Lake,  Minnesota,  July  1987. 

*  (6).  "Ab-Initio  MRD-CI  Calculations  Based  on 
Localized  Orbitals  For  Molecular  Decomposition  and  Reactions  of  Large 
Systems,"  1987  World  Congress  of  World  Association  of  Theoretical  Organic 
Chemistry,  Budapest,  Hungary,  August  1987. 

(7) .  "Ab-Initio  MRD-CI  Calculations  For  Cationic 
Polymerization  of  Oxetanes  Based  on  Localized  Orbitals,"  Division  of 
Physical  Chemistry,  194th  National  Meeting,  American  Chemical  Society,  New 
Orleans,  Louisiana,  August  1987. 

(8) .  "Comparison  of  Ab-Initio  M0DP0T  and  Ab-Initio 
Energy  Partitioned  Potential  Functions  for  Nitromethane  Dimer  Against  Large 
Basis  Set  Calculation,"  Division  of  Physical  Chemistry,  194th  National 
Meeting,  American  Chemical  Society,  New  Orleans,  Louisiana,  August  1987. 


b.  Other  Research  Institutions 


*  (I).  "Ab-Initio  Quantum  Chemical  Calculations  on 
Large  Molecular  Systems  and  Crystals,"  Maryland  Section  American  Chemical 
Society,  Baltimore,  Maryland,  February  1987. 

c.  At  DOD  Meetings  and  Workshops 

*  (1).  "Quantum  Chemical  Characterization  of  Explosive 
Sensitivity,"  Working  Group  Meeting  on  Sensitivity  of  Explosives,  Socorro, 
New  Mexico,  March  1987.  (Presented  by  Dr.  Walter  S.  Koski) 

*  (2).  "Ab-Initio  Quantum  Chemical  Studies  on 
Energetic  Nitrocompounds,"  Sixth  Annual  Working  Group  Meeting  on  Synthesis 
of  High  Energy  Density  Materials,  Concord  Hotel,  Kiamesha  Lake,  New  York, 

May  1987. 

*  (3).  "Ab-Initio  MRD-CI  Calculations  for  Breaking  a 
Chemical  Bond  in  a  Molecule  in  a  Crystal,"  QNR  Workshop  on  Dynamic 
Deformation,  Fracture  and  Transient  Combustion  of  Energetic  Compounds, 

Great  Oak,  Maryland,  May  1987. 

*  (4).  "Desired  Properties  of  Energetic  Compounds  That 
Can  Be  Calculated  Reliably  and  Accurately  From  High  Quality  Ab-Initio 
Quantum  Chemical  Calculations,"  ONR  Workshop  on  Crystalline  and  Polymeric 
Energetic  Materials,  Great  Oak,  Maryland  August  1987. 

2.  To  be  presented 

*  a.  "Ab-Initio  MRD-CI  Calculations  for  Breaking  a  Chemical 
Bond  in  a  Molecule  in  a  Crystal  or  Other  Solid  Environment,"  Sanibel 
International  Symposium  on  Atomic,  Molecular  and  Solid  State  Physics, 
Marineland,  Florida,  March  1988. 

*  b.  "Ab-Initio  MRD-CI  Calculations  for  Breaking  a  Chemical 
Bond  in  a  Molecule  in  a  Crystal  or  Other  Solid  Environment",  Symposium  on 
the  Physics  and  Chemistry  of  Brittle  Fracture,  American  Physical  Society 
National  Meeting,  New  Orleans,  March  1988. 

c.  "Ab-Initio  MRD-CI  Calculations  for  Breaking  a  Chemical 
Bond  in  a  Molecule  in  a  Crystal  or  Other  Solid  Environment,"  3rd  Chemical 
Congress  of  North  America  (joint  with  American  Chemical  Society),  Toronto, 
Canada,  June  1988. 
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3.  Publ ications 

1 .  Already  Publ i shed 

a.  "Ab-Initio  MRD-CI  Calculations  on  the  >C  -  NO^ 

Decomposition  Pathway  of  Nitrobenzene,"  Joyce  J.  Kaufman,  P.  C.  Hariharan, 
S.  Roszak  and  M.  van  Hemert.  An  invited  plenary  lecture  presented  at  the 
Symposium  on  Computational  and  Mathematical  Chemistry,  Can.  Inst.  Chem. 
National  Meeting,  Saskatoon,  Canada,  June  1986.  J.  Comp.  Chem.  8,  736-743 
(1987) . 


b.  "Ab-Initio  Electrostatic  Molecular  Potential  Contour 
Maps  for  Initiation  Step  and  Ab-Initio  MRD-CI  Calculations  for  Propagation 
Step  of  Cationic  Polymerization  of  Oxetanes",  Joyce  J.  Kaufman,  P.  C. 
Hariharan,  S.  Roszak  and  P.  B.  Keegstra.  An  invited  lecture  presented  at 
the  IUPAC  5th  International  Symposium  on  Ring-Opening  Polymerization,  Blois, 
France,  June  1986.  Makromol .  Chem.,  Macromol.  Symposium.  6,  315-330  (1986) 

c.  "Symposium  Note:  More  New  Desirable  Computational 
Strategies  for  Ab-Initio  Calculations  on  Large  Molecules,  Clusters,  Solids, 
and  Crystals",  Joyce  J.  Kaufman,  Int.  J.  Quantum.  Chem.  29,  179-184  (1987) 

d.  "Nonempirical  Atom-Atom  Potentials  for  Main  Components 
of  Intermodular  Interaction  Energy,"  W.  A.  Sokalski,  A.  H.  Lowrey,  S. 
Roszak,  V.  Lewchenko,  J.  M.  Slaisdell,  P.  C.  Hariharan  and  Joyce  J.  Kaufman, 
J.  Comp.  Chem.  7,  693-700  (1986). 

2.  Accepted  for  Publication  and  in  Press 

a.  "Ab-Initio  Potential  Functions  For  Crystals  and  Ab- 
Initio  Crystal  Orbitals,"  Joyce  J.  Kaufman.  An  invited  lecture  presented  at 
the  International  Symposium  on  Molecules  in  Physics,  Chemistry  and  Biology. 
Dedicated  to  Professor  R.  Daudel ,  Paris,  France,  June  1986.  In  press. 
Symposium  Proceedings. 

b.  "Ab-Initio  MRD-CI  Calculations  for  the  Propagation  Step 
in  the  Cationic  Polymerization  of  Oxetanes  Based  on  Localized  Orbitals," 
Joyce  J.  Kaufman,  P.  C.  Hariharan  and  P.  B.  Keegstra.  An  invited  paper 
presented  at  the  International  Sanibel  Symposium  on  Atomic,  Molecular  and 
Solid  State  Theory,  Marineland,  Florida,  March  1987.  In  press,  Int.  J. 
Quantum  Chemistry,  Symposium  Issue. 

c.  "Comparison  of  Ab-Iitio  M0DP0T  and  Ab-Initio  Energy 
Partitioned  Potential  Functions  for  Nitromethane  Dimer  Against  Large  Basis 
Set  Calculations,"  An  invited  paper  presented  at  the  International  Sanibel 
Symposium  on  Atomic,  Molecular  and  Solid  State  Theory,  Marineland,  Florida, 
March  1987.  W.  A.  Sokalski,  P.  C.  Hariharan  and  Joyce  J.  Kaufman.  In 
press,  Int.  J.  Quantum  Chemistry,  Symposium  Issue. 

d.  "Ab-Initio  MRD-CI  Calculations  on  Protonated  Cyclic 
Ethers.  I.  Protonation  Pathways  Involve  Multipotential  Surfaces 
(Protonation  of  Oxetane)  II.  Differences  from  SCF  in  Dominant 
Configurations  Upon  Opening  Non-Protonated  Oxirane  Rings  (epoxides),"  An 
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invited  lecture  presented  at  the  International  Sanibel  Symposium  on  Atomic, 
Molecular  and  Solid  State  Theory,  Marineland,  Florida,  March  1987.  Joyce  J 
Kaufman,  P.  C.  Hariharan,  S.  Roszak  and  P.  B.  Keegstra,  In  press,  Int.  J. 
Quantum  Chemistry,  Symposium  Issue. 
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VI  Project  Personnel 


Joyce  J.  Kaufman,  Ph.D. 


Principal  Investigator 


P.  C.  Hariharan,  Ph.D. 


Research  Scientist 


Overall  responsibility  for  implementing  new  program  developments  and 
conversion  to  Cray  computers.  Quantum  chemical  calculations  on  energetic 
polymers,  MRO-CI,  GAMESS  and  POLY-CRYST  calculations 


Philip  B.  Keegstra,  Ph.D. 


Postdoctoral  (October  1986  -  August  1987) 


MRD-CI  calculations  on  cationic  polymerization  of  oxetanes, 
implementation  of  inclusion  of  multipoles  into  POLY-CRYST,  assistance  on 
implementing  multipoles  into  MRD-CI  and  with  breaking  chemical  bond 
calculations,  test  calculations  on  POLY-CRYST 


S.  Roszak,  Ph.D. 


Visiting  Scientist  (February  1987  -  November  1987) 


Implementation  of  multipoles  into  MRD-CI  and  carrying  out  such 
calculations  for  breaking  a  chemical  bond  in  a  molecule  in  a  crystal,  MRD-CI 
calculations  on  propagation  step  of  cationic  polymerization  of  oxetanes. 


W.  A.  Sokalski,  Ph.D. 


Visiting  Scientist  (July  -  August  1987) 


Inclusion  of  correlation  in  calculation  of  multipoles  and  use  in 
intermolecular  calculations. 
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